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SUMMARY 


Measurements  of  acoustic  signals  in  the  Boston  area  during 
summer  1979  are  reported.  The  measurements  demonstrate  that  the 
secondary  sonic  boom  events  frequently  reported  by  New  England 
residents  during  the  summer  are  created  by  Concorde  flights  off 
the  New  England  coast  en  route  to  Kennedy  Airport  in  New  York  City. 
Reception  times  for  such  events  are  highly  correlated  with  Concorde 
flight  times;  the  time  delay  between  aircraft  time  at  closest  point 
of  approach  to  Hyannis  MA  and  signal  reception  at  Malden  MA  is 
constant  within  53  seconds.  Ray  trace  computations,  based  on  best 
available  meteorological  data  and  on  Federal  Aviation  Administration 
(FAA)  radar  measurements  of  Concorde  flight  profiles,  predict  signal 
arrival  times  within  20  seconds  of  that  observed  and  predict  wave 
front  azimuthal  angles  within  4°. 

Signal  amplitudes  show  wide  fluctuations  from  f 1 ight - to - f 1 ight 
and  from  day-to-day.  The  average  measured  event  maximum  peak-to- 
peak  pressure  amplitude  for  May,  June,  July,  August,  and  the  first 
half  of  September  are  0.15,  0.14,  0.19,  0.20,  and  0.11  pounds  per 
square  foot  (lbs/ft2).  The  largest  individual  values  are  0.63, 

0.69,  0.65,  0.59,  and  0.35  for  the  five  months,  respectively. 
Simultaneous  measurements  carried  out  at  other  locations  in  the 
Boston  area  show  that  there  are  3lso  wide  variations  with  geographic 
location  for  the  same  event.  A  brief  set  of  measurements  for  four 
days  in  July  were  made  at  App lebachsv i 1 le  PA.  These  ’Measurements 
show  similar  day-to-day  variability  and  are  correlated  with  Concorde 
flights  into  Dulles  Airport. 

The  ray  trace  computations  indicate  that  the  events  with 
larger  amplitudes  are  associated  with  strong  focusing  of  the 
secondary  boom  at  the  reception  site.  Although  the  locations  where 
strongest  focusing  occur  vary  significantly  with  meteorological 
conditions,  the  ray  theory  computations  show  that  they  generally 
fall  within  the  greater  Boston  area  during  the  summer  months. 


xii 


The  data  taken  during  this  measurement  program  and  the  theoret 
ical  interpretation  of  the  results  provide  the  basis  for  evaluating 
the  effectiveness  of  alternate  operational  procedures  to  reduce  the 
frequency  and  intensity  of  secondary  sonic  booms  within  the  U.S. 
mainland . 


1.  INTRODUCTION 


This  report  documents  the  results  of  a  Secondary  Sonic  Boom 
Detection  and  Assessment  program  conducted  in  New  England  during  the 
summer  of  1979  by  the  U.S.  Department  of  Transportation  -  Trans¬ 
portation  Systems  Center  (DOT/TSC)  for  the  Office  of  Environment 
and  Energy  of  the  Federal  Aviation  Administration. 

Secondary  sonic  boom  events  were  reported  in  the  New  England 
area  during  the  summer  of  1978.  These  events  were  reported  as 
muffled  "thumps"  and  low-frequency  "rumbles"  and  occurred  pre¬ 
dominantly  in  the  morning  hours  between  8:00  and  10:30  am.  Reports 
were  received  from  citizens  as  far  north  as  Maine  and  along  New 
England  coastal  regions  of  New  Hampshire,  Massachusetts,  and  Rhode 
Island.  Because  of  the  widespread  nature  of  the  reported  events  and 
their  low  frequency  nature,  it  was  speculated  at  that  time  that 
these  events  originated  from  the  shock  wave  of  supersonic  aircraft 
operating  off  the  New  England  coast. 

Preliminary  investigation  of  the  phenomena  by  DOT/TSC  during 
September  1978  showed  that  there  was  some  correlation  in  timing 
between  boom  data  measured  north  of  Boston  and  visual  sightings  of 
the  New  York-bound  Concorde  at  Chatham  MA.  However,  at  that  time 
data  was  insufficient  to  clearly  establish  the  Concorde  as  the 
primary  source  of  these  secondary  sonic  booms. 

A  measurement  program  was  thus  devised  to  determine  conclusively 
the  source  of  the  energy  that  precipitated  the  public  reports  and  to 
provide  an  overall  assessment  of  the  propagation  and  disturbance 
mechanism . 


2.  TECHNICAL  BACKGROUND 


FAA  flight  rules  require  supersonic  transports  (SSTs)  to  fly  at 
subsonic  speeds  over  U.S.  land  areas  in  order  to  prevent  sonic  booms 
from  impacting  the  U.S.  environment.  For  aircraft  approaching  U.S. 
boundaries,  flight  rules  specify  operational  procedures  be  such  that 
direct  sonic  shock  waves  do  not  encroach  upon  the  U.S. 

Although  these  requirements  are  adequate  to  prevent  primary 
shock  waves  from  impacting  U.S.  population  centers,  secondary  long 
distance  effects  are  knovm  to  occur.  ’  *  A  principal  propagation 
mechanism  causing  such  long  distance  effects  is  refraction  caused 
by  wind  and  temperature  gradient  effects  at  altitudes  between  20  and 
60  km  (the  mesosphere).  Sound  rays  that  carry  upward  traveling  sonic 
booms  to  such  altitudes  can  be  bent  back  toward  the  ground  if  these 
gradients  cause  the  effective  sound  speed  to  increase  with  height; 
the  basic  mechanism  has  been  understood  for  some  time  in  connection 
with  the  abnormal  zones^  of  audibility  associated  with  the  sound 
heard  at  large  distances  from  gunfire  and  explosions.  Such  down¬ 
ward  refraction  can  also  take  place  in  the  t hermosphere^ ’ ^  (altitudes 
above  80  km);  but  the  high  attentuation  and  the  lengthening  of  the 
shock  duration  at  such  high  altitudes  would  render  such  thermospheric 
refracted  arrivals  much  less  likely  to  be  audible  by  the  time  they 
return  to  the  ground. 

The  secondary  booms  that  arrive  via  a  refracted  path  that  bends 
downward  in  the  mesosphere  can  be  of  two  basic  types  (see  Figure  1). 
Type  I  originates  in  the  shock  wave  that  is  generated  by  the  air¬ 
plane  and  initially  travels  obliquely  upward.  Type  II  is  initiated 
from  a  downward  shock  wave  that  reflects  from  the  sea  surface 
relatively  close  to  the  aircraft's  trajectory,  then  travels  up  to 
the  thermosphere  and  returns  to  the  ground.  In  either  case,  it  is 
possible  for  the  propagation  to  be  such  that  perceivable  acoustic 
events  occur  in  restricted  on-shore  regions.  This  refractive  mode 
of  low  frequency  acoustic  propagation  extends  the  possible  impact 
region  of  an  SST's  supersonically  generated  acoustical  disturbance 


beyond  its  primary  boom  carpet.  The  ranges  at  which  such  secondary 
paths  may  reach  the  ground  can  be  anywhere  from  50  to  200  miles 
beyond  the  primary  carpet. 


Meteorological  conditions  play  an  important  role  in  this  long- 
range  propagation  of  acoustical  shock  waves.  Of  the  family  of 
rays  that  leave  the  aircraft  and  are  refracted,  only  a  smrll 
fraction  reach  the  ground.  Only  when  the  sum  of  the  component  of 
the  wind  velocity  in  the  direction  of  propagation  and  the  sound 
speed  (at  some  altitude  above  the  flight  altitude)  exceeds  the 
corresponding  sum  at  the  ground,  is  it  possible  for  these  rays  to 
reach  the  ground  beyond  the  primary  carpet  area.  The  geometry  of 
the  paths  may  be  such  that  large  areas  on  the  ground  are  not 
reached  by  any  rays  (creating  quiet  zones)  ,  while  in  areas  further 
from  the  flight  trajectory,  rays  may  touch  down,  such  that  energy 
is  received.  This  phenomenon,  known  as  anomalous  or  abnormal  sound, 
has  been  studied  on  various  occasions  in  connection  with  large 
explosions  near  the  ground.  Up  to  now  it  has  not  been  studied  in 
much  detail  in  connection  with  sonic  boom  propagation.  The 
extreme  sensitivity  of  the  propagation  paths  to  the  sound  speed 
and  wind  velocity  profiles  and  large  scale  turbulence  in  the 
atmosphere  suggests  that  the  energy  received  at  any  point  will 
vary  substantially  with  meteorological  conditions. 

This  effect  of  upper  atmospheric  wind  currents  explains  the 
absence  of  reports  in  New  England  during  the  winter  months  when  the 
stratospheric  winds  are  predominantly  westerly,  thus  inhibiting 
refractive  propagation.  On  the  other  hand,  during  the  summer 
months  when  strong  easterly  winds  predominate,  the  secondary  sonic 
boom  events  become  evident. 

An  important  consquence  of  these  meteorological  influences 
is  that  the  energy  reaching  the  ground  via  such  long  paths  is  un¬ 
equally  spread  over  the  ground  and  that  it  is  virtually  impossible 
to  predict  precisely  where  the  energy  will  be  most  concentrated. 
Given  the  atmosphere's  basic  structure,  only  a  relatively  small 
fraction  of  the  acoustic  energy  created  by  the  Concorde’s  super¬ 
sonic  flight  should  reach  the  ground  via  paths  that  are  refracted 


downwards  in  the  upper  atmosphere.  If  this  energy  were  uniformly 
spread  over  a  wide  range  of  lateral  distances,  then  it  is  doubtful 
that  the  secondary  booms  would  ever  be  audible.  However,  focusing 
can  produce  an  audible  effect  in  a  relatively  small  area.  Thus, 
the  reported  events  in  New  England  are  caused  by  circumstances  when 
wind  and  temperature  profiles,  aircraft  heading,  speed,  and  loca¬ 
tion  in  the  right  combination  produce  an  augmentation  (or  focusing) 
of  acoustic  energy  at  the  location  of  the  observer. 


3,  EXPERIMENTAL  APPROACH 


3.1  BASIC  APPROACH 


Low-frequency  pressure  changes  were  measured  with  three  micro¬ 
phone  measuring  systems  deployed  in  a  triangular  array  during  the 
period  May  3,  l':l79  to  September  14  ,  1979.  The  signals  from  each 
of  the  microphones  along  with  a  time  code  signal  were  recorded  on 
a  multitrack,  FM  tape  recorder.  The  differences  in  arrival  time 
among  microphones  and  the  position  and  dimensions  of  the  array 
were  used  to  compute  the  direction  of  the  arriving  wave  front. 

This  arriving  ray  vector  was  then  used  to  project  back  to  the 
signal  source. 


In  addition  to  the  three  microphone  array  deployed,  a  single 
microphone  recording  system  with  time  code  synchronization  was  also 
deployed  on  occasion  at  several  sites  throughout  the  Greater  Boston 
area.  Measurements  were  made  simultaneously  with  the  multi¬ 
microphone  and  the  single  microphone  systems  to  show  the  vari¬ 
ability  of  the  signals  measured  with  site  location. 


3.2  MEASUREMENT  SITES 

A  prime  measurement  site  was  selected  within  the  expected 
"first  bounce"  impact  area  in  the  Middlesex  Fells  Reservation  in 
Malden  MA  (Figure  2) .  A  lightly  wooded  area  was  selected  to  pro¬ 
vide  natural  wind  reduction  with  relatively  little  attenuation  of 
the  expected  signals. 

Six  other  sites  were  chosen  within  a  25  mile  radius  of  Boston 
to  show  the  variability  of  the  data  with  site  location.  These  sites, 
Wilmington,  Georgetown,  Marlboro,  Medfield,  Sharon,  and  Cohasett 
MA,  are  shown  in  Figure  2. 

A  measurement  site  was  also  sought  with  coordinates  relative 
to  the  flight  track  of  inbound  Goncorde  flights  to  Dulles  Inter¬ 
national  Airport  similar  to  those  of  the  Malden  site  relative  to 


inbound  Concorde  flight  to  Kennedy  Airport  (JFK) .  Measurements 
were  made  during  the  week  of  July  23,  1979  at  a  site  in  Applebachs- 
ville  PA  that  met  this  requirement  (see  Figure  3) . 


3.3  INSTRUMENTATION  AND  DEPLOYMENT 

Low-frequency  (infrasonic)  microphone  systems,  consisting  of 
B$K  4146  one-inch  condenser  microphones  and  BSK  2631  Microphone 
Carrier  System  were  used  (Figure  4).  The  infrasonic  data  was  band- 
filtered  (0.5  to  30  Hz)  and  recorded  on  a  four-track  FM  Instrumenta 
tion  Tape  Recorder  HP  Model  3960A.  The  microphones  were  enclosed 
in  a  weather  proof  wooden  enclosure  (30"  x  20"  x  16")  to  further 
reduce  the  effects  of  wind  noise.  An  acoustic  calibration  signal 
of  125  Hertz  (Hz)  at  a  level  of  114  decibels  (db)  re  20  micro- 
Pascal  was  recorded  on  tape  to  provide  a  reference  level  for  the 
data  reduction  instrumentation. 

Three  low-frequency  microphone  systems  were  placed  in  a  tri¬ 
angular  array  250  feet  apart  (Figure  4).  This  triad  was  deployed 
at  the  Malden  MA  site  except  during  the  weeks  of  July  9-13,  1979 
and  July  23-26,  1979,  when  it  was  deployed  at  Georgetown  MA  and 
Applebachsville  PA  respectively. 

A  single  low-frequency  microphone  system  was  deployed  at  the 
Malden  MA  site  during  these  two  weeks  to  provide  continuity  of  the 
data.  Occasionally,  the  single  system  was  also  deployed  at  sites 
in  Wilmington,  Marlboro,  Medfield,  Sharon,  and  Cohasett  MA  to  pro¬ 
vide  an  indication  of  the  variability  of  the  secondary  boom  that 
occurs  within  a  25  mile  radius  of  Boston. 

A  vibration-measuring  system  consisting  of  an  Endevco  Acceler¬ 
ometer  Model  2217E  with  signal  conditioners  was  deployed  in  a  res¬ 
idential  home  in  Wilmington  MA  on  June  20,  1979.  The  accelerometer 
was  cemented  to  the  center  of  a  54  by  60  inch  by  one-quarter  inch 
plate  glass  window  located  in  a  room,  facing  east,  at  the  rear  of 
the  house.  Lateral  acceleration  levels  (2  Hz  to  1  kilohertz  (kHz)) 
of  the  glass  were  recorded  in  synchronism  with  the  infrasonic  data 
from  the  low-frequency  microphone  system  located  outside  approxi¬ 
mately  100  feet  from  the  window. 
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To  measure  and  record  events  as  perceived  by  a  listener,  an 
acoustic  measuring  system  using  a  half  inch,  General  Radio  1962- 
9610  Electret  condenser  microphone  was  also  deployed  on  occasion. 
Measured  outdoor  acoustic  data  in  the  frequency  band  15  Hz  to 
20  kHz  was  recorded  on  a  Nagra  model  IVSJ  instrumentation  tape 
recorder . 

A  time  codi  signal  was  recorded  on  all  recorders  for  exact  time 
synchronization. 

3.4  FLIGHT  DATA 

Flight  track  data  was  obtained  from  the  FAA  Data  Systems  Office 
in  Nashua  NH  from  the  air  traffic  control  radar  system.  This  data 
includes  aircraft  flight  number,  speed,  altitude,  position,  and  time 
of  day.  Physical  data  from  the  on-board  recorders  of  several 
selected  flights  was  obtained  from  British  Airways  and  Air  France 
to  supplement  the  radar  data.  (See  typical  tracks  to  JFK  and  Dulles 
Airports  in  Figures  2  and  3,  respectively.) 

In  addition,  an  aircraft  receiver  tuned  to  Boston  Air  Route 
Traffic  Control  Center  was  used  to  assist  in  aircraft  identification. 

3.5  WEATHER  DATA 

Local  weather  was  recorded  continuously  on  site.  Niigh  altitude 
data  were  obtained  from  rocketsonde  observations  from  Wallops 
Island  VA.  Local  high  altitude  data  were  obtained  from  balloon 
observations  and  from  satellite-infrared  radiance  data  -  all  obtained 
through  the  courtesy  of  the  Air  Resources  Laboratory  of  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA) . 

3.6  DATA  ANALYSIS 

Utilizing  a  Nicolet  Scientific  Dual  Channel  Fast  Fourier 
Transform  (FFT)  Analyzer  Model  411A,  cross-correlation  techniques 
yielded  time  differences  and  chronology  for  the  signals  received 
by  the  microphones  of  the  triangular  array.  Since  the  microphones 
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were  deployed  in  a  triangular  array  and  one  side  of  the  triangle 
ran  north  to  south,  this  data  could  be  used  to  calculate  the  azimuth 
and  elevation  angles  of  the  arriving  ray  vector  as  follows: 

-1  At21 

Azimuth  angle  (ft)  =  tan  -r-z — 

At  23 

where  At2i  =  time  difference  between  microphone  2  and  1 
At 2 3  =  t ime  difference  between  microphone  2  and  3 


Elevation  angle  =  cos 


-1  cAt23 
d  cosO 


where  c  =  local  speed  of  sound 

d  =  the  distance  between  microphones  2  and  .3 

The  narrow-band  frequency  spectra  and  peak  pressure  of  the 
signals  were  also  obtained.  The  absolute  level  of  this  data  was 
referenced  to  the  recorded  calibration  signal,  i.e.,  114  dB  sound 
pressure  level  corresponds  to  a  pressure  of  0.33  peak  lbs/ft^. 

The  narrow  band  frequency  spectrum  was  referenced  to  a  narrow 
band  measurement  (1.25  Hz  line  width)  of  the  recorded  calibra¬ 
tion  signal  (114  dB  re  20  microPascal). 

Graphic  history  recordings  of  the  measured  infrasonic  neak 
pressure  changes  were  made  using  a  SAM- EL  Instruments  Recti  - 
graph  model  8S.  (Selected  representative  events  recorded  are  in¬ 
cluded  in  Section  4.2.) 


4.  MEASUREMENT  DATA 


During  the  period  May  3,  1979  to  September  14,  1979,  infra- 
sonic  measurements  were  made  at  Malden  MA,  at  six  other  sites  in  the 
Greater  Boston  area  (See  Figure  2) ,  and  at  a  site  in  Applebachs- 
ville  PA  (see  Figure  3).  Measurements  were  made  on  a  daily  basis 
throughout  the  period  at  the  Malden  site  with  measurements  simul¬ 
taneously  made  (one  site  at  a  time)  at  each  of  the  other  sites  to 
obtain  a  measure  of  the  variability  of  the  data  with  geographic 
location. 

4 . 1  SUMMARY  DATA 

Tables  1  through  ID  contain  summary  data  of  the  events  re¬ 
corded  on  a  daily  basis  at  the  Malden  site  for  May,  June,  July, 
August,  and  September  1979,  respectively.  Table  IE  summarizes  the 
data  measured  at  the  remaining  seven  sites. 

An  inspection  of  these  tables  shows  the  variability  of  the 
measured  peak-to-peak  pressure  change  resulting  from  the  unstable 
stratospheric  winds.  Of  the  events  measured,  seven  were  recorded 
with  peak-to-peak  pressure  changes  in  excess  of  0.5  lbs/ft2.  The 
largest  signal  (0.69  lbs/ft2)  was  measured  on  June  14,  1979.  The 
lowest  signals  were  measured  during  September,  1979,  and  signaled 
the  end  of  the  1979  season  of  easterly  st ratospheric  winds  and 
the  uninhibited  propagation  of  secondary  sonic  booms  that  were 
generated  off  shore.  The  1979  measurement  program  was  therefore 
terminated  on  September  14,  1979. 

On  September  6,  1979,  Air  France  Flight  AF-001  flew  a  track 
approximately  40  miles  south  of  the  usual  JFK-bound  flight  track. 

The  new  flight  track,  operational  procedures,  and  the  absence  of 
easterly  stratospheric  winds  on  that  particular  day  resulted  in 
the  low  pressure  levels  measured  for  that  particular  flight. 

A  comparison  of  data  in  Table  IE  (Miscellaneous  Sites)  with 
data  at  the  Malden  site  during  the  same  time  frame  shows  the  vari¬ 
ability  of  the  data  because  of  geographic  location.  It  is  seen 
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that  the  levels  measured  in  Georgetown  were  generally  greater  than 
those  in  Malden.  At  the  remaining  five  sites  in  Massachusetts, 
the  reverse  was  generally  true.  This  is  further  illustrated  in 
Section  4.2  in  the  comparison  of  pressure  time  histories. 

Angular  dimensions  of  the  arriving  ray  vector  are  also  in¬ 
cluded  in  these  tables.  The  azimuth  angle  is  the  calculated  heading 
of  the  arriving  ray  vector  measured  clockwise  from  true  north.  The 
measured  azimuth  angle  averaged  approximately  278  degrees  re  true 
north  at  Malden.  If  the  aircraft  is  taken  to  have  a  250  degree 
heading,  the  azimuth  angle  indicates  (see  Section  4.4,  Ray  Tracing) 
the  aircraft  was  traveling  at  approximately  Mach  1.31  when  the  shock 
wave  was  launched  (see  Figure  2)  . 

The  elapsed  time  tabulated  is  the  difference  between  the  time 
of  the  measured  maximum  peak-to-peak  pressure  change  and  the  time 
of  the  closest  point  of  approach  (CPA)  of  the  aircraft  to  Hyannis 
MA,  as  determined  by  radar.  The  statistical  consistency  of  the 
elapsed  time  over  the  summer  indicates  conclusively  that  these  JFK 
bound  flights  are  the  source  of  the  secondary  shock  waves  in  New 
England.  Further,  it  indicates  that  the  shock  waves  In  question 
are  launched  from  a  specific  portion  of  the  flight  track  of  the 
aircraft . 

Technicians  on  site  heard  the  majority  of  the  signals  measured 
and  in  some  cases  were  able  to  obtain  a  sound  level  meter  reading 
(on  the  fast  scale)  unencumbered  by  local  disturbances.  These  have 
also  been  tabulated.  The  highest  sound  pressure  level  (SPL) 
measured  on  June  18,  1976,  overloaded  a  sound  level  meter  set  for 
a  full  scale  reading  of  60  dBA.  A  peak-to-peak  pressure  change  of 
0.68  lbs/ft2  was  recorded  for  this  event,  one  of  the  highest  levels 
measured  during  the  summer. 

4.2  PRESSURE  SIGNATURES 

Representative  time  histories  of  pressure  changes  recorded  are 
presented.  Figures  5  to  16  show  the  day-to-day  as  well  as  the 
hour-to-hour  (or  f 1 ight-to- f 1 ight )  variability  of  the  signals 
measured  at  the  Malden  measurement  site  during  the  1979  measurement 
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program.  Extremes  in  the  hour-to-hour  variability  can  be  seen  in 
Figures  10  and  12  for  July  18  and  26,  1979,  while  Figure  13  for 
August  15,  1979  shows  a  marked  similarity  in  the  two  events 
measured . 

An  examination  of  the  Malden  signatures  shows  that  not  only  do 
the  amplitudes  change  but  the  complexity,  duration,  and  multiplicity 
of  the  waveforms  also  vary.  The  waveforms  appear  less  complex  in 
May  and  September  than  those  in  the  intervening  months.  Azimuth 
calculations  of  various  portions  of  the  pressure  signatures  show 
the  waveforms  are  made  up  of  signals  arriving  from  a  variety  of 
azimuth  angles.  When  the  upper  level  winds  begin  to  inhibit  the 
refracted  wave  from  touching  down,  as  in  early  May  and  late 
September,  the  number  of  paths  are  reduced  and  the  pressure 
signatures  become  less  complex  (see  Figures  15  and  16). 

A  direct  comparison  can  be  made  in  Figures  17  through  24  of 
the  signals  measured  at  Malden  with  the  simultaneous  measurements 
recorded  at  Marlboro,  Medfield,  Cohasset,  Georgetown,  and 
Sharon  MA.  Both  the  amplitude  and  shape  of  the  signatures  are  seen 
to  be  significantly  different.  The  larger  of  the  levels  was 
measured  in  Malden,  with  the  exception  of  those  measured  in 
Marlboro  and  Georgetown. 

Figures  25  and  26  show  pressure  time  histories  of  data  recorded 
in  Applebachsville  PA  from  Concorde  flights  in-bound  to  Dulles  Inter¬ 
national  Airport.  The  signatures  shown  have  the  same  degree  of  day- 
to-day  variability  as  have  the  Malden  data.  The  phenomenon  is  similar 
to  that  observed  in  New  England. 

Figure  27  contains  time  histories  of  both  exterior  infrasonic 
pressure  changes  and  window  vibration  accelerations  measured  on 
June  20,  1979,  at  a  residential  home  in  Wilmington  MA.  The  lateral 
accelerations  of  the  54  by  60  by  one-quarter  inch  plate  glass 
window  are  seen  to  be  in  synchronism  with  the  outside  pressure 
changes  resulting  from  the  secondary  sonic  boom.  In  this  instance, 
technicians  both  inside  and  outside  the  home  heard  the  event. 


Figure  28  presents  expanded  pressure  time  histories  of  two 
reasonably  uncomplicated  sonic  signals  measured  on  May  3  and  14, 

1979.  In  addition,  synchronous  histories  are  shown  in  Figure  28A 
of  the  pressure  measured  with  the  triangular  array  of  infrasonic 
sensors.  Note  that  the  differences  between  times  of  occurrence 
are  clearly  discernable  in  the  three  time  histories. 

The  time  difference  between  the  traces  for  microphone  2  and 
microphone  1  shows  that  the  wave  front  arrived  from  an  easterly 
direction,  i.e.,  the  wave  front  arrived  at  microphone  1  prior  to 
microphone  2.  In  addition,  the  signal  arrived  at  microphone  2 
prior  to  microphone  3  (as  shown)  ,  which  indicates  a  ray  direction 
from  the  south.  This  combined  information  shows  the  signal  arrived 
roughly  from  a  southeasterly  direction.  Actual  geometric  calcula¬ 
tions  of  time  differences  (obtained  more  precisely  by  cross -correlat ion 
techniques)  shows  che  azimuth  angle  or  heading  of  the  arriving  ray 
vector  to  be  in  this  case  279  degrees  corrected  to  true  north. 

Pressure  time  histories  of  selected  more  complex  signals  are 
shown  in  Figures  29-33.  Each  figure  supplements  the  overall  pressure 
time  history  of  the  event  with  expanded  segments  of  the  history, 
showing  the  fine  detail  of  the  complex  pressure  changes. 

4.3  FREQUENCY  SPECTRA 

A  narrow-band  frequency  analysis  was  performed,  averaging  a 
4-second  period  of  data  centered  at  the  time  of  the  maximum  recorded 
pressure  change.  The  frequency  spectra  for  several  representative 
events  are  presented  in  Figures  34-41.  The  signal  drop-off  above 
30  Hz  is  characteristic  of  the  TSC  infrasonic  measuring  system, 
which  was  adjusted  for  a  high  frequency  cutoff  at  30  Hz.  The 
narrow  band  frequency  spectrum  of  the  acoustic  event  recorded  on 
May  14,  1979  is  presented  in  Figure  35  for  comparison  with  the 
narrow  band  spectrum  of  the  infrasonic  signal  measured  simultaneously 
(Figure  34B) .  Note  that  Figure  35A  shows  the  presence  of  acoustic 
energy  up  to  approximately  100  Hz.  An  expansion  of  the  low- 
frequency  portion  of  this  spectra  is  shown  in  Figure  3SB.  It  is 
seen  to  be  an  extension  of  the  data  measured  by  the  infrasonic 
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system  (Figure  34B) .  Note  that  the  low-frequency  roll-off  of  the 
data  in  Figure  35B  is  characteristic  of  the  acoustic  measuring 
system,  which  has  a  flat  response  above  20  Hz  and  is  approximately 
6  dB  down  at  15  Hz.  The  measured  "A"  weighted  SPL  of  this  signal 
was  46  dBA  in  an  environment  with  low  level  ambient  of  40  dBA. 

Narrow  band  analysis  of  various  other  acoustic  signals  measured 
during  the  program  (acoustic  instrumentation  was  not  deployed  every 
day)  did  not  exhibit  any  significant  increase  in  the  high  frequency 
(greater  than  100  Hz)  portion  of  the  spectra. 

Figures  36-39  contain  the  narrow  band  frequency  spectra  for 
representative  events  measured  in  Malden  MA  during  June,  July, 
August,  and  September,  1979.  (See  Figures  8,  9,  13  and  16  for 
pressure  time  history  data.)  Figure  40  contains  frequency 
data  for  events  measured  on  July  11,  1979,  in  Georgetown  MA  and  can 
be  compared  with  the  frequency  spectra  (Figure  37)  for  the  same 
events  measured  in  Malden  MA.  (Also  see  the  pressure  histories  in 
Figures  22  and  9.) 

Figure  41  contains  frequency  spectra  for  signals  from  two 
flights  measured  in  Applebachsvi lie  PA.  (See  Figures  25  and  26 
for  the  pressure  time  histories). 

Figure  42  contains  the  frequency  spectra  for  the  infrasonic 
event  measured  on  June  20,  1979  in  Wilmington  MA,  and  the  spectra 
of  the  window  pane  vibrations.  (See  Figure  27  for  the  corresponding 
time  histories.) 
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5.  PUBLIC  REPORTS  AND  ACOUSTIC  DATA 


Several  hundred  reports  were  received  by  the  FAA  from  citizens 
in  the  Greater  Boston  area  of  sounds  indicative  of  secondary  sonic 
booms  during  the  summer  of  1978.  The  Planning  and  Appraisal  Staff 
of  the  New  England  Region  FAA  called  several  of  these  citizens  dur¬ 
ing  the  latter  part  of  May,  1979,  to  solicit  their  assistance  in 
observing  the  phenomena  and  to  keep  records  of  their  observations 
for  correlation  with  the  measurements  program.  Table  2  includes 
observations  from  five  citizens  plus  some  unsolicited  reports 
received  during  May,  1979,  by  the  FAA  New  England  Region. 

Note  that  the  reported  observation  times  coincide  very  closely 
with  the  measured  data  as  shown  in  Table  2.  Reports  received  by 
the  FAA  from  May  through  July,  1979,  constitute  the  basis  for  Figure 
43,  that  shows  geographically  the  area  from  which  the  reports  were 
received. 

Citizen  reports  indicate  the  secondary  booms  were  "heard." 

TSC  techncians  at  the  measurement  site  indicate  they  "heard"  many 

of  the  secondary  booms  recorded.  As  shewn  in  the  Summary  Data 

(Tables  1-1E),  the  technicians  "heard"  the  secondary  booms  when 

the  maximum  peak-to-peak  pressure  changes  exceed  approximately 
2 

0.1  lb/ft  .  On  several  days,  sound  level  meter  measurements 
were  made  on  site  and  levels  up  to  58  dBA  (fast  scale)  were 
measured.  On  one  occasion  (July  18,  1970)  the  sound  level  meter 
overloaded  and  indicated  a  level  in  excess  of  60  dBA.  It  should 
be  noted  from  Table  IB  that  one  of  the  largest  peak-to-peak  pres¬ 
sures  measured  was  recorded  that  day  (0.68  lbs/ft^). 

Frequency  analysis  of  these  data  in  narrow  band  (see  Figure 
35)  and  in  one  third  octave  (data  not  shown)  shows  little,  if  any, 
energy  above  a  level  of  30  to  40  dB  re  20  micro  Pascal  (uPa)  in 
the  frequency  range  above  100  Hz.  One  third  octave  analysis  was 
done  using  a  Gen  Rad  One  Third  Octave  Real  Time  Analyzer. 
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As  shown  in  the  expanded  pressure  time  history,  Figure  33D, 
the  rise  time  of  the  signal  is  less  than  10  milliseconds.  For  an 
impulse,  this  would  indicate  a  cutoff  frequency  of  approximately 
100  Hz,  which  agrees  with  the  above  frequency  analysis. 
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6.  DATA  INTERPRETATION 


Many  of  the  features  of  the  data  measured  in  the  field  exper¬ 
iments  can  be  interpreted  using  basic  acoustical  principles. 

6.1  WAVEFORM  SEGMENTS 

While  the  pressure  time  histories  recorded  in  the  field  show 
a  great  amount  of  variability,  it  is  possible  with  guidance  from 
theory  and  ray  computations  to  perceive  three  different  types  of 
arrivals  in  most  of  the  time  histories.  The  data  (Figures  8B  and 
31)  received  at  Malden  that  was  caused  by  the  BA-171  flight  (see 
Figure  44)  of  June  20,  provides  an  illustration  of  such  a  classifica¬ 
tion.  The  segments  B  and  C  in  Figure  31A  correspond  to  a  ground 
wave  and  to  a  type  I  secondary  sonic  boom ,  respectively.  The 
segments  D,  E,  and  F  correspond  collectively  to  a  type  II  secondary 
sonic  boom. 

The  ground  wave  is  of  very  small  amplitude  and,  in  many  of  the 
studied  events,  is  either  too  small  to  be  detected  or  else  requires 
careful  scrutiny  of  the  recording  before  it  is  found.  It  is  always 
the  first  arrival  and  typically  arrives  about  1  minute  before  the 
type  I  secondary  boom.  The  propagation  path  of  the  ground  wave  goes 
from  the  flight  trajectory  of  the  airplane  to  the  edge  of  the 
primary  carpet  along  a  direct  (but  refracted)  geometrical  acoustic 
ray.  However,  at  the  edge  of  the  primary  carpet,  the  energy 
carried  by  impinging  direct  rays  splits.  Most  of  the  energy  re¬ 
flects  back  into  the  atmosphere,  but  a  portion  continues  to  propa¬ 
gate  along  the  ground  in  the  same  general  direction  as  a  creeping 
wave.  The  creeping  wave  continuously  sheds  energy  into  the 
atmosphere  and  therefore  dies  out  exponentially  with  increasing 
propagation  distance.  Consequently,  very  little  amplitude  remains 
by  the  time  it  reaches  the  Malden  site.  The  theory  of  this  ground 

7 

wave  dates  back  to  fundamental  papers  by  C.  L.  Pekeris  (1946)  and 

g 

by  D.C.  Pridmore-Brown  and  U.  Ingard  (1955).  The  application  of 
the  theory  to  sonic  booms  has  been  subsequently  discussed  by 

9 

R.  Onyeonwu  (1975). 
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The  type  I  secondary  sonic  boom  arrival  travels  along  a  geo¬ 
metric  acoustic  ray  path  (see  Figures  1  and  44)  that  initially  goes 
obliquely  upward  from  the  flight  trajectory,  but  bends  downward  as 
it  propagates  upward  until  eventually  the  ray  is  turned  back  toward 
the  ground.  Typically  the  turning  point  is  at  45,000  to  55,000 
meters  in  altitude.  Amplitudes  at  Malden  of  type  I  arrivals  are 
highly  variable,  but  usually  the  type  I  arrival  is  weaker  than  the 
type  II  arrival.  The  tentative  explanation  for  this  tendency  is 
that  the  energy  associated  with  shocks  propagating  obliquely  down¬ 
wards  from  a  Concorde  in  supersonic  flight  is  substantially  more 
than  that  associated  with  shocks  propagating  obliquely  upwards. 

The  type  II  secondary  sonic  boom  waveform  segment  typically 
arrives  about  30  to  50  seconds  after  the  type  I  segment.  The 
delay  is  because  it  travels  a  longer  path,  going  from  aircraft 
trajectory  to  ground  (impacting  within  primary  carpet)  ,  then  re¬ 
flecting  and  propagating  up  to  the  stratosphere  and  mesosphere, 
and  refracting  downward  back  to  the  ground. 

6.2  EFFECT  OF  UPPER  ATMOSPHERIC  WINDS 

The  simplest  model  of  sonic  boom  propagation  is  based  on 
geometrical  acoustics.  This  model  predicts  type  I  and  type  II 
sonic  boom  arrivals  beyond  the  primary  carpet  area  only  when  the 
sum  of  the  wind  velocity  component  in  the  direction  of  propaga¬ 
tion  and  the  sound  speed  at  an  altitude  above  the  flight 
altitude,  exceeds  the  corresponding  sum  at  the  ground.  (See 
Figure  45.)  This  implies  that  the  presence  in  the  Boston  area  of 
secondary  booms  caused  by  inbound  Concorde  flights  is  strongly 
affected  by  the  east-west  component  of  the  wind  velocity  at  heights 
of  30  to  60  kilometers  (km) .  Strong  winds  blowing  east  to  west  at 
such  altitudes  should  be  associated  with  the  presence  of  secondary 
booms.  The  very  small  amplitudes  received  during  early  May  can 
be  explained  by  observing  meteorological  conditions  of  the  upper 
atmospheric  winds  that  are  blowing  in  the  wrong  direction  or  are 
very  light  at  the  time. 
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To  test  the  hypothesis  just  mentioned  the  scatter  plot  in 
Figure  46  was  prepared.  Each  point  corresponds  to  an  event  recorded 
at  the  Malden  site  during  May  or  June.  The  amplitudes  that  con¬ 
stitute  the  vertical  coordinates  are  the  peak-to-peak  maxima  listed 
in  Tables  1  and  1A.  The  horizon  coordinates  are  east-to-west  wind 
components  at  55,000  meter  altitude  for  the  corresponding  day  of 
flight.  The  latter  numbers  come  from  estimates  based  on  meteoro¬ 
logical  data  obtained  from  satellites;  these  estimates  are  for  the 
Boston  area  and  are  supplied  to  us  by  the  NOAA  Air  Resources 
Laboratory  (see  Appendix  A) .  The  accuracy  of  such  meteorological 
estimates  is  not  known,  so  the  scatter  plot  also  constitutes  a 
rough  test  of  the  veracity  of  the  NOAA  estimates. 

Figure  46  demonstrates  that  the  probabiliiy  of  receiving  a 
large  amplitude  event  is  small  unless  the  east-to-west  wind  speed 
at  the  considered  altitude  exceeds  16  meters  per  second  (m./s). 

There  is  one  large  amplitude  event  (associated  with  the  BA-171 
flight  of  May  9)  that  is  a  drastic  exception  to  this  rule,  but 
the  day  to  day  reliability  of  the  NOAA  predictions  is  unknown. 

It  is  expected,  however,  that  statistically,  they  are  valid 
predict  ions . 

Even  when  the  upper  atmospheric  winds  are  of  high  amplitude 
and  blowing  in  the  right  direction,  there  is  no  binding  theoretical 
reason  that  the  amplitudes  of  secondary  booms  at  a  specified 
measurement  site  should  be  large.  The  detailed  geographical  dis¬ 
tribution  of  the  received  secondary  boom  energy  has  a  complicated 
and  even  somewhat  erratic  dependence  on  the  details  of  the 
meteorological  structure  all  along  the  propagation  path  as  well 
as  on  the  flight  profile,  so  any  simple  rules  for  predicting 
large  amplitude  secondary  sonic  booms  must  be  regarded  in  a 
statistical  sense.  From  this  point  of  view,  Figure  46  would 
appear  to  give  substantial  confirmation  of  the  hypothesis  that 
occurrence  of  such  events  is  strongly  associated  with  stronger 
magnitude  east-to-west  winds  in  the  upper  atmosphere. 


6.3  MEASURED  AND  PREDICTED  ARRIVAL  TIMES 


Given  sufficient  meteorological  data  and  sufficient  data  con¬ 
cerning  the  aircraft's  flight  profile,  it  is  possible  to  compute 
predictions  for  the  arrival  times  of  the  ground  wave  and  of  the 
type  I  and  type  II  secondary  sonic  boom  waveforms  at  any  given 
measurement  site.  The  meteorological  data  needed  consists  of 
temperature,  wind  speed,  and  wind  direction  versus  height  up  to 
an  altitude  of  at  least  50  km.  Some  sample  profiles  are  shown  in 
Figures  47,  48,  and  49.  Flight  profile  data  required  consists  of 
aircraft  position,  altitude,  speed,  and  heading  versus  time. 

Meteorological  profiles  for  the  months  of  May,  June,  and  July 
were  furnished  to  TSC  by  the  NOAA  Air  Resources  Laboratory;  flight 
profiles  (radar  trackings  of  inbound  Concordes)  were  furnished  by 
the  FAA  (see  Appendix  A  for  a  fuller  discussion).  From  the  subset 
of  events  for  which  both  meteorological  data  and  flight  profile 
data  were  available,  the  thirteen  listed  in  Table  3  were  selected 
for  detailed  computational  study. 

With  the  exception  of  two  subtleties  (described  below),  the 
calculation  of  arrival  times  at  the  Malden  site  involves  a  straight- 
foward  (although  nontrivial)  application  of  sonic  boom  and  geo¬ 
metrical  acoustics  principles.  Rays  are  continously  being  "shed" 
along  the  flight  profile  and  they  obey  certain  constraints  in 
their  possible  initial  directions  that  are  related  to  the  air¬ 
craft's  heading  and  Mach  number.  Subsequent  propagation  of  the 
shedded  rays  is  affected  by  the  meteorological  profiles  and  fol¬ 
lows  the  basic  rules  of  geometrical  acoustics.  The  computer 
program  written  by  TSC  for  this  purpose  is  summarized  in  Appendix  A. 

One  of  the  problems  that  had  to  be  overcome  was  that  the  Malden 
site  was  invariably  too  close  to  the  flight  profile  for  the  geo¬ 
metrical  acoustics  model  per  se  to  predict  the  arrival  time  of 
the  type  II  secondary  sonic  boom.  All  of  the  computed  ray  paths 
that  left  the  flight  profile  would  fall  beyond  Malden.  The  miss 
distance  of  Malden  from  the  nearest  such  type  II  ray  impact  point, 
however,  was  typically  not  large,  being  usually  only  about  10 
to  20  percent  of  the  total  propagation  distance. 


The  geometrical  acoustics  model  of  sonic  boom  propagation 
is  approximate  and  valid  at  best  only  in  the  limit  of  very  high 
frequencies  and  negligible  attenuation,  so  it  is  not  surprising 
to  find  a  type  II  arrival  at  Malden  when  the  simplistic  model 
doesn't  yield  any  ray  path  connecting  the  flight  trajectory  with 
Malden.  A  proper  theoretical  interpretation  of  the  geometrical 
acoustics  theory  is  that  it  furnishes  the  basic  "super-structure" 
of  the  wavefield  to  which  the  "billowing  sails”  of  the  full-wave 
theory  are  "sewn." 

It  was  observed  from  the  computed  results  that,  at  those  more 
distant  points,  where  the  geometrical  acoustics  model  did  predict 
both  type  I  and  type  II  ray  arrivals,  the  time  difference  in  their 
arrivals  along  a  given  azimuth  direction  was  nearly  independent 
of  range.  Consequently,  it  was  argued  that  the  same  difference 
would  apply  to  the  Malden  site.  The  computed  type  II  arrival 
times  at  Malden  shown  in  Table  3  are  therefore  sums  of  the  cor¬ 
responding  computed  type  I  arrival  times  and  increments  based  on 
parameters  of  type  II  rays  that  reach  further  points. 

Another  modification  made  in  the  computation  was  to  assume 
that  winds  had  no  effect  on  the  ground  wave  transit  speed,  so  the 
appropriate  speed  is  simply  the  sound  speed  on  the  ground.  The 
reason  for  leaving  out  the  wind  speed  is  that  the  NOAA  estimates 
of  wind  velocities  are  based  on  the  geostrophic  assumption  rather 
than  direct  measurements.  Such  an  assumption,  although  good  at 
moderate  to  high  altitudes,  is  inappropriate  near  the  ground  and 
within  the  earth's  boundary  layer.  With  regard  for  the  frictional 
effect  of  the  ground  and  in  absence  of  other  data,  zero  wind  speed 
at  the  ground  seems  to  be  the  most  rational  assumption  of  comparable 
s impl icity . 

The  measured  arrival  times  with  which  the  computed  arrival 
times  in  Table  3  are  compared  are  times  of  maximum  amplitude  for 
the  corresponding  waveform  segment.  In  some  cases,  the  actual 
waveform  portion  corresponding  to  a  given  propagation  mode  had  a 
duration  of  as  much  as  30  seconds.  If  a  single  time  must  be 
selected  to  compare  with  theory,  however,  the  time  of  maximum 
amplitude  should  be  least  biased. 
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The  discrepancies  between  computed  and  measured  times  are 
very  small  compared  to  a  representative  travel  time,  aircraft 
trajectory  to  Malden,  of  an  arriving  acoustic  ray.  The  latter 
is  of  the  order  of  750  to  800  seconds.  In  contrast,  the  average 
(absolute)  discrepancies  in  Table  3  are  12  seconds,  18  seconds,  and 
24  seconds  for  the  ground  wave,  type  I  and  type  II  arrivals,  re¬ 
spectively.  (The  latter  are  1.6  percent,  2.3  percent,  and  3.0 
percent  of  800  seconds.) 

The  agreement  strongly  supports  identification  of  Concorde 
flights  bound  for  JFK  as  the  cause  of  the  events  recorded  at 
Malden.  It  also  supports  the  suppositions  that  divide  the  total 
waveforms  into  ground  wave,  type  I,  and  type  II  segments. 

The  arrival  time  agreements  do  not,  however,  constitute  a 
sensitive  test  of  the  accuracy  of  the  meteorological  data  because 
the  travel  time,  equal  to  path  integral  of  reciprocal  of  wave  speed, 
tends  to  be  insensitive  to  minor  fluctuations  in  atmospheric 
properties.  Neither  do  they  provide  a  thorough  test  of  the  ac¬ 
curacy  of  ray  path  predictions  because  Fermat's  principle 
guarantees  that  acoustic  travel  time  along  a  path  connecting  two 
points  should  be  insensitive  to  small  variations  in  the  path. 

The  May  30  meteorological  data  was  such  that  no  type  I  or 
type  II  ray  paths  are  predicted  that  touch  the  ground  anywhere  to 
the  northwest  of  the  flight  track.  This  does  not  necessarily  mean 
that  such  arrivals  are  precluded  by  the  full  wave  theory,  but  the 
development  of  a  scheme  for  estimating  their  arrival  times  did  not 
seem  warranted  in  the  present  study.  The  small  amplitudes  re¬ 
corded  for  the  AF  001  event  for  this  particular  date  suggest, 
however,  that  the  meteorological  conditions  were  not  favorable 
for  the  occurrence  of  strong  secondary  sonic  booms. 


6.4  RAY  AZIMUTHS  AND  ELEVATIONS 

The  measured  azimuths  and  elevations  (shown  in  Tables  1-1E) 
of  the  arriving  rays  correspond  to  the  waveform  segment  with 
greatest  amplitude.  Usually,  as  explained  in  Section  6.1,  this 
is  the  type  II  arrival  waveform  segment.  The  arrival  azimuths  for 


successive  waveform  segments  were  measure!  for  a  number  of  events; 
the  results  show  some  variation  hut  the  type  I  arrival  character¬ 
istically  arrives  with  an  azimuth  2°  less  than  that  of  the  type  II 
arrival.  (Smaller  azimuths  correspond  to  arrivals  that  leave  the 
flight  track  earlier.) 

Our  computer  predictions  for  points  where  both  types  of  rays 
are  expected  according  to  the  geometrical  acoustics  model  are  con¬ 
sistent  with  this  trend;  differences  for  points  just  to  the  west 
of  Malden  are  typically  1.5°  to  2°. 

The  plot  in  Figure  50  of  computed  versus  measured  arrival  ray 
azimuth  angles  takes  the  computed  numbers  to  be  those  for  the  type  I 
rays,  since  the  geometrical  acoustics  model  usually  does  not  yield 
such  a  number  for  type  II  arrivals  at  the  Malden  site.  Adding  2° 
to  each  computed  azimuth  (as  might  be  suggested  by  the  remarks  in 
the  preceding  paragraph)  would  not  change  the  slope  of  the  linear 
least  squares  fit  line,  but  would  raise  it.  This  would  change 
the  point  of  its  intersection  to  282°  instead  of  278°.  From 
either  point  of  view,  the  agreement  seems  substantial. 

Table  4  lists  computed  and  measured  elevation  angles  for  the 
events  selected  for  analysis.  Here  again,  the  computed  numbers 
correspond  to  type  I  rays  arriving  at  Malden;  the  measured  values 
are  based  for  the  most  part  on  the  local  horizontal  transit  speed 
of  the  type  II  waveform  segment  across  the  array  at  the  Malden 
site.  The  measured  ray  elevation  angle  is  the  arc-cosine  of  the 
ratio  of  the  sound  speed  at  the  ground  to  the  measured  transit 
speed  across  the  array.  The  computed  elevation  angles  are  derived 
from  the  geometrical  acoustics  ray  paths  computed  by  the  program. 

Because  the  elevation  angles  are  small  and  because  the  cosine 
(approximately  1)  of  a  small  angle  is  relatively  insensitive  to 
small  variations,  the  measured  elevation  angles  are  probably  not 
too  accurate.  Also,  the  computed  elevation  angles  are  very 
sensitive  to  small  variations  in  the  assumed  meteorological 
model's  sound  speed  and  wind  velocity  on  the  ground.  For  these 
reasons,  the  poor  correlation  of  measured  and  computed  values  in 
Table  4  is  to  be  expected.  The  results,  however,  are  of 
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comparable  magnitudes  and  consequently  support  the  hypothesis  that 
the  secondary  booms  at  the  Malden  site  have  traveled  along  paths 
that  reached  stratospheric  altitudes  before  returning  to  the  ground. 

A  related  comparison  of  measured  and  computed  quantities  is 
shown  in  Figure  51;  local  horizontal  transit  speeds  ( phase  ve¬ 
locities)  of  waveforms  across  the  microphone  array  are  predicted 
by  the  basic  theory  to  be  the  same  as  the  sum  (at  the  highest 
altitude  reached  by  the  incident  ray)  of  the  sound  speed  and  wind 
velocity  component  in  propagation  direction.  This  prediction  is  in 
accord  with  the  extension  of  Snell's  law  that  takes  winds  into 
account.  The  comparison  of  the  measured  and  computed  local  hori¬ 
zontal  transit  speeds  should  therefore  be  insensitive  to  the  local 
meteorological  conditions  and  topography  at  the  array  site;  the 
comparison  tests  instead  the  reliability  of  the  meteorological  pro¬ 
files  at  turning  point  altitudes  (35  km  to  50  km  above  the  earth's 
surface) . 

6.5  SECONDARY  BOOM  FOCUS  LINES 

The  geometrical  acoustics  theory  predicts,  for  type  I  secondary 
boom  arrivals,  a  line  ( secondary  boom  focus  line)  that  divides  re¬ 
gions  where  no  type  I  rays  hit  the  ground  from  regions  where  they 
do  hit  the  ground.  An  analogous  line  is  predicted  that  separates 
regions  struck  by  type  II  geometrical  acoustics  rays  from  regions 
where  they  do  not  strike.  (The  reasons  for  referring  to  such  lines 
as  focus  lines  are  stated  further  below.) 

Figure  44  shows  the  two  secondary  boom  focus  lines  computed 
by  the  TSC  computer  program  (see  Appendix  A)  for  the  BA  171  flight 
of  June  30.  Similar  computational  results  are  depicted  in 
Figures  52,  53,  54,  and  55  for  the  May  31  BA  171  flight,  the 
June  27  AF  001  flight,  and  the  July  18  AF  001  and  BA  171  flights. 

The  type  I  secondary  boom  focus  line  is  always  closer  to  the 
flight  track  than  is  the  type  II  focus  line.  The  separation 
distance  between  the  two  lines  is  typically  about  30  nautical 
miles.  The  locations  and  shapes  of  the  two  lines  depend  on  the 
details  of  the  aircraft's  flight  profile  as  well  as  on  the 
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meteorological  data.  That  small  variations  in  flight  profile  max- 
affect  the  lines,  even  when  the  meteorology  does  not  change,  is 
evident  from  a  comparison  of  Figures  54  and  55. 

For  all  the  events  analyzed,  the  predicted  type  II  focus  line 
always  lies  further  from  the  flight  track  than  does  the  Malden 
site.  The  Malden  site  typically  falls  between  the  two  lines,  but 
in  some  instance^  the  type  I  focus  line  passes  through  some  point 
west  of  Malden.  1 igure  54  for  the  July  18  AF-001  flight  exempli¬ 
fies  the  latter  possibility.  Far  to  the  east-north-east  of  the 
Boston  area,  the  two  lines,  as  indicated  in  Figure  44,  eventually 
approach  asymptotes  that  are  parallel  to  the  flight  track  and 
correspond  to  the  region  of  level  rectilinear  flight  at  a  constant 
Mach  number  of  2. 

It  is  possible  that  type  I  and  type  II  lines  exist  south  of 
the  flight  track  but  the  meteorological  conditions  (upper  atmosphere 
winds  blowing  east  to  west)  for  the  considered  events  are  likely 
to  prevent  the  secondary  sonic  boom  ray  from  returning  to  the 
ground  south  of  the  flight  track.  Our  computations,  however,  have 
been  confined  to  disturbances  in  the  New  England  area. 

The  southernmost  tips  of  the  type  I  and  type  II  focus  lines 
that  appear  in  the  figures  are  associated  with  the  last  point  in 
the  corresponding  FAA  radar  track  tabulation  for  which  the  air¬ 
craft  was  flying  at  supersonic  speed  relative  to  the  ground.  The 
FAA  tabulations  are  at  time  intervals  of  12.5  seconds,  so  any 
extensions  of  the  secondary  boom  focus  lines  to  points  further 
south  must  correspond  to  less  than  12.5  seconds  of  flight  (or  to 
less  than  2  miles  of  flight).  The  lines,  if  extended,  would  resemble 
arcs  of  circles  with  radii  of  approximately  200  km  and  260  km, 
respectively,  extending  down  to  the  flight  track.  However,  because 
the  secondary  boom  energy  generated  per  unit  time  is  roughly  con¬ 
stant  and  because  the  energy  generated  during  the  last  few  seconds 
is  spread  over  such  a  wide  geographical  area,  the  amplitudes  at 
points  south  of  where  the  computed  lines  stop  are  expected  to  be 
relatively  low. 
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The  identification  of  the  computed  lines  as  focus  lines  results 
because  they  describe  lines  along  which  the  geometrical  acoustics 
rays  converge.  The  ray  density  (number  of  rays  per  unit  area) 
formally  become  infinite  and  the  ray  tube  area  goes  to  zero  along 
such  lines.  The  geometrical  acoustics  theory  breaks  down  at  these 
focus  lines  but  predicts  large  amplitudes  near  the  lines.  On  the 
illuminated  side  the  amplitude  should  fluctuate  somewhat  because 
of  interference  between  rays  that  have  touched  and  have  not  yet 
touched  the  caustic  surface,  but  the  general  trend  is  for  amplitude 
to  decrease  as  the  inverse  fourth  root  of  distance  from  the  focus 
line . 

On  the  side  of  the  focus  line  where  the  geometrical  acoustics 
theory  predicts  an  absence  of  impinging  rays,  the  full  wave  theory^ 
predicts  that  the  amplitude  dies  off  with  distance  y  from  the  focus 
line  roughly  as 

exp  [-(4TT/3)(2/X2R)1/2(tane)3/2  y3/2] 

where  X  is  the  wavelength,  R  is  a  distance  of  the  order  of  200  km, 
and  the  angle  0  is  of  the  order  of  10°.  The  distance 

L  =  (3/4Tr)2/3(A2R/2)1/3/tane 

can  be  taken  as  a  measure  of  the  half-width  of  the  strip  over  which 
the  influence  of  the  focus  extends  on  either  side.  For  X  =  340  m 
(f  =  1  Hz),  L  is  of  the  order  of  5  km.  (The  choice  of  1  Hz  for 
a  representative  frequency  follows  from  an  examination  of  the 
expanded  pressure  time  histories  in  Figures  28-33.) 

The  hypothesis  that  larger  amplitudes  are  associated  with  the 
proximity  of  the  measurement  site  to  a  focus  line  is  tested  by 
the  scatter  plots  in  Figures  56  and  57. 

The  former,  Figure  56,  plots  the  measured  maximum,  peak-to- 
peak,  pressure  change  for  the  type  I  waveform  segments  versus 
computed  perpendicular  distance  of  the  Malden  site  from  the  type  I 
focus  line  (which  is  the  focus  line  closer  to  the  flight  track). 
Here,  positive  distance  implies  Malden  lies  on  the  far  side  of 
the  focus  line  in  relation  to  the  flight  track  or,  equivalently, 
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that  t ho  focus  lino  lies  between  Malden  and  flu*  flight  track, 
for  simplicity,  t ho  least  squares  fit  line,  which  is  included  to 
indicate  the  general  trend  rather  than  to  approximate  the  functional 
dependence  predicted  hv  theory,  is  based  on  the  assumption  that 
maximum  amplitude  is  symmetrical  about  the  focus  line. 

The  general  trend  supports  the  hypothesis,  but  there  is  con¬ 
siderable  scatter.  Our  interpretation  of  this  scatter  is  that  the 
amplitude  dependence  on  distance  from  a  focus  line  must  be  viewed 
in  terms  of  probabilities.  Proximity  to  the  focus  line  increases 
the  probability  that  a  high  amplitude  will  be  measured,  but  does 
not  guarantee  it.  One  reason  for  the  absence  of  a  deterministic 
relationship  is  because  focusing  is  also  affected  by  turbulence** 
and  by  departures  of  the  atmosphere  from  perfect  stratification. 

Another  source  of  fluctuations  is  that  our  predictions  of 
focus  line  locations  are  of  dubious  accuracy.  The  meteorological 
information  available  is  not  sufficiently  detailed  to  allow  such 
lines  to  be  pinpointed.  A  rough  guess  is  that  the  accuracy  in 
focus  line  prediction  is  of  the  order  of  4  20  km.  If  this  is  so, 
then  the  actual  sharpness  of  the  focusing  will  be  considerably 
blurred  when  one  plots  amplitude  versus  computed  distance  to  focus 
line  rather  than  versus  true  distance  to  focus  line.  From  this 
perspective,  the  rather  weak  slope  of  the  linear  least  squares 
fit  line  in  Figure  56  is  not  surorising,  even  though  the  character¬ 
istic  scale  for  the  width  of  the  focus  region  that  is  predicted  by 
theory  is  only  of  the  order  of  ’0  km. 

The  preceding  remarks  apply  also  for  the  plot  in  Figure  57, 
in  which  maximum  peak-to-peak  pressure  changes  for  type  II 
waveform  segments  arc  plotted  versus  distance  to  the  type  II 
secondary  boom  focus  line.  Here,  however,  the  positive  sense  of 
direction  is  opposite  to  that  adopted  in  Figure  56.  For  all 
the  events  studied,  Malden  was  always  closer  to  the  flight  track 
than  was  the  computed  type  fl  focus  lino.  A  positive  distance  of, 
for  example,  50  km  in  Figure  57  means  the  type  II  focus  line  passed 
through  a  point  northwest  of  Malden  that  was  roughly  30  km  distant 
from  the  measurement  site.  Also,  because  the  type  II  focus  line 


20 


is  more  distant  from  the  track  than  is  the  type  I  focus  line,  it 
is  expected  that  its  computed  location  is  even  less  precise.  The 
estimate  is  that  the  accuracy  in  focus  line  prediction  in  this 
case  is  of  the  order  of  *  30  km. 
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7.  CONCLUSIONS 


Low  frequency  acoustic  signals  arrived  at  regular  intervals  in 
the  New  England  area  during  the  late  spring  and  summer  months  of 
1979.  Arrival  times  and  detailed  analysis  of  the  measured  signals 
leave  no  doubt  that  they  are  caused  by  inbound  Concorde  flights 
into  JFK  International  Airport  in  New  York  City.  During  every 
period  when  the  measurement  system  was  in  operation  from  May  3 
to  September  14,  1979,  and  when  a  Concorde  was  in  an  inbound 
approach  path  to  JFK,  a  signal  was  recorded.  A  similar  statement 
applies  to  the  measurements  at  Applebachsville  PA  (July  23  to 
27,  1979)  in  regard  to  Concorde  flights  into  Dulles. 

For  the  New  England  events,  the  measured  arrival  times  of  peak 
signal  amplitudes  correlated  extremely  well  with  time  of  Concorde 
passage  to  the  closest  point  of  approach  to  Hyannis  MA.  Knowledge 
of  the  latter  gives  the  former  with  an  rms  error  of  only  53 
seconds.  Furthermore,  if  account  is  taken  of  the  detailed  flight 
profile  of  the  aircraft,  the  event  arrival  times  can  be  predicted 
to  within  20  seconds  using  acoustic  propagation  theory  and  available 
meteorological  data. 

During  May,  55  percent  of  the  measured  events  at  Malden  were 
audible  to  the  field  technicians.  In  June,  94  percent  were  audible; 
in  July,  9"  percent;  in  August,  86  percent;  in  the  first  half  of 
September,  53  percent.  Consequently,  it  is  concluded  that,  for 
the  approach  profiles  currently  being  flown,  almost  every  Concorde 
flight  into  JFK  during  the  summer  months  generates  audible  sound 
in  the  Boston  area.  The  four  measurements  at  Applebachsville 
PA  during  July  suggest  a  similar  conclusion  concerning  flights 
into  Dulles. 

The  preponderance  of  reports  during  the  summer  months  is  well 
explained  in  terms  of  acoustic  propagation  theory  by  the  seasonal 
changes  in  stratospheric  wind  direction.  During  summer,  the  winds 
in  the  altitude  region,  30  km  to  60  km,  blow  east  to  west  and  cause 
secondary  sonic  booms  propagating  high  above  the  aircraft  flight 
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altitude  to  refract  back  down  to  the  ground.  Detailed  calculations 
based  on  radar  trackings  of  Concordes  and  on  best  available 
meteorological  data  for  13  events  yield  respectable  agreement 
with  data  in  regard  to  arrival  times,  azimuth  angles,  and  local 
horizontal  transit  speeds.  Therefore,  the  basic  mechanisms  by 
which  the  secondary  booms  are  conveyed  to  the  Boston  area  can  be 
regarded  as  well  understood. 

The  computations  and  the  data  yield  the  conclusions  that  the 
signals  received  in  the  Boston  area  from  the  Concorde  flights  arc 
generated  during  the  last  minute  of  supersonic  flight,  when  the 
aircraft  is  off  the  New  F.ngland  coast  between  68°  and  69°  IV  longitude. 

Although  there  are  noticeable  similarities  among  the  waveforms 
received  from  different  flights  and  on  different  dates,  there  are 
wide  variations  in  details  of  signature  shape  and  in  amplitudes. 

The  maximum  (within  a  given  waveform)  pressure  change,  peak-to-peak , 
for  events  measured  at  Malden  in  May  ranged  from  0.02  to  0.63 

lbs/ft2;  in  July,  they  ranged  from  0.04  to  0.68  lbs/ft";  in 

?  ‘  * 

August,  they  ranged  from  0.07  to  0.59  lbs/ft";  in  September, 
they  ranged  from  less  than  0.01  to  0.32  lbs/ft2. 

The  large  amplitude  events  are  ascribed  to  an  uneven  distribu¬ 
tion  of  acoustic  energy  impinging  on  the  ground;  a  region  receiving 
an  abnormally  large  amplitude  is  a  region  of  acoustic  focusing. 

Just  where  such  focusing  occurs  is  highly  sensitive  to  details  in 
the  atmosphere's  meteorological  structure,  but  the  basic  theoretical 
model  that  neglects  longitudinal  and  latitudinal  variations  of 
weather  profiles  predicts  two  focus  lines  along  which  abnormally 
high  amplitudes  may  be  expected.  Both  lines  tvpically  pass  through 
the  Boston  area  with  a  separation  distance  of  about  30  nautical 
miles.  The  Malden  measurement  site  typically  lies  betAveen  the 
focus  lines,  but  small  variations  in  flight  profile  and  in 
meteorological  structure  can  cause  the  focus  lines  to  move  in 
either  direction  by  as  much  as  35  nautical  miles.  The  comparison  of 
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measured  peak  amplitudes  w  i  t  h  distance's  of  the  Malden  site  from  the 
computed  locus  lines  supports  the  hypothesis  that  larger  amplitudes 
are  caused  bv  focusing. 

The  data  taken  during  t h i s  measurement  program  and  the 
theoretical  i  nterpret  a  t  ion  of  the'  results  [)rovide  the  basis  for 
evaluating  the  effectiveness  <  f  alternate  operational  procedures 
to  reduce  the  frequency  and  intensity  of  secondary  sonic  booms 
within  the  U.S.  mainland. 


FIGURE  2.  FLIGHT  TRACK  -  JFK-BOUND  FLIGHT  BRITISH  AIRWAYS  FLIGHT  BA-171  JUNE  20,  1979 


FIGURE  4.  SYSTEM  SCHEMATIC  DIAGRAM 
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FIGURE  5.  PRESSURE  TIME  HISTORIES  SOURCE:  BRITISH  AIRWAYS 
FLIGHT  BA- 171  MALDEN  MA 
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FIGURE  7.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  JUNE  13,  1979 
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FIGURE  8.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  JUNE  20,  1979 


PRESSURE  CHANGE  -  (pounds  per  square  foot) 


A)  Source:  Air  France  Flight  AF-001 

B)  Source:  British  Airways  Flight  BA-171 


FIGURE  9.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  JULY  11,  1979 
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PRESSURE  CHANGE  -  (pounds  per  square 


FIGURE  10.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  JULY  18,  1979 


44 


PRESSURE  CHANGE  -  (pounds  per  square  foot) 
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FIGURE  11.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  JULY  23,  1979 
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PRESSURE  CHANGE  -  (pounds  per  square  foot) 
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FIGURE  12.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  JULY  26,  1979 
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A)  Source:  Air  France  Flight  AF-001 

B)  Source:  British  Airways  Flight  BA-171 


FIGURE  13.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  AUGUST  15,  1979 
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A)  Source:  Air  France  Flight  AF-001 

B)  Source:  British  Airways  Flight  BA-171 


FIGURE  14.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  AUGUST  21,  1979 
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FIGURE  16.  PRESSURE  TIME  HISTORIES  MALDEN  MA  -  SEPTEMBER  12,  1979 
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FIGURE  17.  PRESSURE  TIME  HISTORIES  SOURCE:  BRITISH  AIRWAYS 
FLIGHT  NO.  BA-171  JUNE  27,  1979 
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FIGURE  18.  PRESSURE  TIME  HISTORIES  SOURCE:  BRITISH  AIRWAYS 
FLIGHT  NO.  BA-171  JUNE  28,  1979 


52 


PRESSURE  CHANGE  -  (pounds  per  square  foot) 
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FIGURE  19.  PRESSURE  TIME  HISTORIES  SOURCE:  BRITISH  AIRWAYS 
FLIGHT  NO.  BA-171  JUNE  29,  1979 
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PRESSURE  CHANGE  -  (pounds  per  square  foot) 
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FIGURE  20.  PRESSURE  TIME  HISTORIES  SOURCE:  AIR  FRANCE  FLIGHT 
NO.  AF-001  JULY  9,  1979 
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FIGURE  21.  PRESSURE  TIME  HISTORIES  SOURCE:  BRITISH  AIRWAYS 
FT.IGHT  NO.  BA-171  JULY  9,  1.979 
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PRESSURE  CHANGE  -  (pounds  per  square  foot) 
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FIGURE  22.  PRESSURE  TIME  HISTORIES  GEORGETOWN  MA 
(See  Figure  9  for  July  11,  1979  data  at  Malden  MA) 
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PRESSURE  CHANGE  -  (pounds  per  square 
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FIGURE  23.  PRESSURE  TIME  HISTORIES  SOURCE:  AIR  FRANCE  FLIGHT 
NO.  AF-001  AUGUST  2,  1979 
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FIGURE  24.  PRESSURE  TIME  HISTORIES  SOURCE:  BRITISH  AIRWAYS 
FLIGHT  NO.  BA-171  AUGUST  2,  1979 
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FIGURE  25.  PRESSURE  TIME  HISTORIES  APPLEBACHSVILLE  PA 
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FIGURE  27.  TIME  HISTORIES  WILMINGTON  MA  -  JUNE  20,  1979 
SOURCE:  BRITISH  AIRWAYS  FLIGHT  NO.  BA-171 
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FIGURE  29.  EXPANDED  PRESSURE  TIME  HISTORY  MALDEN  MA  -  JUNE  13 
1979  SOURCE:  AIR  FRANCE  FLIGHT  AF-001  (CONT.) 
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FIGURE  30.  EXPANDED  PRESSURE  TIME  HISTORY  MALDEN  MA  -  JUNE  20 
1979  SOURCE:  AIR  FRANCE  FLIGHT  AF-001 
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FIGURE  31.  EXPANDED  PRESSURE  TIME  HISTORY  MALDEN,  MA  -  JUNE  20, 
1979  SOURCE:  BRITISH  AIRWAYS  FLIGHT  BA-171 
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FIGURE  31.  EXPANDED  PRESSURE  TIME  HISTORY  MALDEN  MA  -  JUNE  20, 
1979  SOURCE:  BRITISH  AIRWAYS  FLIGHT  BA-171  fCONT . ) 
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PRESSURE  CHANGE  -  (pounds  per  square  foot) 
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FIGURE  32.  EXPANDED  PRESSURE  TIME  HISTORY  MALDEN,  MA  -  JUL 
1979  SOURCE:  BRITISH  AIRWAYS  FLIGHT  BA-171  (CONT.) 
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PRESSURE  CHANGE  -  (pounds  per  square  foot) 


A)  Overall  Time  History 


B-E)  Expanded  Time-History  Segments 

B)  0939:51.8  -  0939:57.8 

C)  0940:41.2  -  0940:47.2 

D)  0940:59.6  -  0941:05.6 

E)  0941:11.6  -  0941:17.6 


FIGURE  33.  EXPANDED  PRESSURE  TIME  HISTORY  MALDEN,  MA  -  JULY  18, 
1979  SOURCE:  BRITISH  AIRWAYS  FLIGHT  BA-171 
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PRESSURE  CHANGE  -  (pounds  per  square  foot) 
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FIGURE  33.  EXPANDED  PRESSURE  TIME  HISTORY  MALDEN,  MA  -  JULY  18, 
1979  SOURCE:  BRITISH  AIRWAYS  FLIGHT  BA-171  (CONT.) 
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RMS  SOUND  PRESSURE  LEVEL  -dB  re  20  uPa 


A)  Source:  British  Airways  Flight  BA-171 

May  3,  1979 

B)  Source:  British  Airways  Flight  BA-171 

May  14,  1979 


FIGURE  34.  INFRASONIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT, 
MAI. DEN  MA 

(See  Figure  28  for  Pressure  Time  Histories) 
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A)  Source:  British  Airways  Flight  BA-171 

B)  Source:  British  Airways  Flight  BA-171 


FIGURE  35.  ACOUSTIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT,  MALDEN  MA 
MAY  14,  1979 


RMS  SOUND  PRESSURE  LEVEL  -  dB  re  20  uPa 
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A)  Source:  Air  France  Flight  AF-001 

B)  Source:  British  Aiiwrays  Flight  BA-171 


FIGURE  36.  INFRASONIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT,  MALDEN  MA 
JUNE  20,  1979 

(See  Figure  30  and  31  for  Pressure  Time  Histories) 
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A)  Source:  Air  France  Flight  AF-001 

B)  Source:  British  Airways  BA-171 


FIGURE  38.  INFRAS0NIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT,  MALDEN  MA 
AUGUST  15,  1979 

(See  Figure  13  for  Pressure  Time  Histories) 


(0 

cu 

a 

o 

fN 


(1) 

a 

'O 


a 

w 

> 

w 

p 

w 

« 

p 

m 

in 

w 

pi 

cu 

P 


o 

in 


in 


A)  Source:  Air  France  Flight  AF-001 

B)  Source:  British  Airways  Flight  BA-171 


FIGURE  39.  INFRASONIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT,  MAIDEN  MA 
SEPTEMBER  12,  1979 

(See  Figure  16  for  Pressure  Time  Histories'! 
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A)  Source:  Air  France  Flight  AF-001 

B)  Source:  British  Airways  Flight  BA-171 


FIGURE  40.  INFRASONIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT,  GEORGETOWN 
MA  JULY  11,  1979 

CSee  Figure  22  for  Pressure  Time  Histories) 
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A)  Source:  Air  France  Flight  AF-053 

July  23,  1979 

B)  Source:  British  Airways  Flight  BA-189 

July  26,  1979 


FIGURE  41.  INFRASONIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT, 
APPLEBACHSVI LLE  PA 

(See  Figures  25  and  26  for  Pressure  Time  Histories 
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A)  Infrasonic  Exterior  pressure  Change 

B)  Window  Vibrational  Acceleration 


FIGURE  42.  INFRASONIC  FREQUENCY  SPECTRA  ACOUSTIC  EVENT,  WILMINGTON 
MA  SOURCE:  BRITISH  AIRWAYS  FLIGHT  BA-171  JUNE  20,  1979 

(See  Figure  27  for  Level  Time  Histories) 
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FIGURE  44.  SECONDARY  SONIC  BOOM  FOCUS-LINES  SOURCE:  BRITISH  AIRWAYS  FLIGHT 
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FIGURE  45.  EFFECTIVE  SOUND  SPEED  PROFILE  BOSTON  MA 
JUNE  20,  1979 

CSee  Figures  47,  48,  and  49  for  meteorological  profiles) 


(•)  Rocketsonde  Data-Wallops  Island,  VA . 

(x)  Rawi nsonde (Balloon  borne)  Data-Boston,  MA . 
(o)  Satellite  Radiance  Data-Boston,  MA . 


FIGURE  47.  TEMPERATURE  PROFILE  JUNE  20,  1979 
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FIGURE  50.  COMPUTED  VERSUS  MEASURED  ARRIVAL  RAY  AZIMUTH  ANGLES  MALDEN  MA 
(42°26'34"N,  71°05’07"W) 

(See  Table  4  for  event  a-m  Identification) 
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FIGURE  51.  COMPUTED  VERSUS  MEASURED  HORIZONTAL  TRANSIT  SPEED  MALDEN  MA 


FIGURE  52.  SECONDARY  SONIC  BOOM  FOCUS-LINES  SOURCE:  BRITISH  AIRWAYS  FLIGHT  BA-171 
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FIGURE  56.  MAX  PRESSURE  CHANGE  VERSUS  DISTANCE  MALDEN  MA  TO  TYPE  1  SECONDARY  BOOM 

FOCUS-LINE  . 

(See  Table  3  for  event  a-m  identification) 


TABLE  1 

SECONDARY  SONIC  BOOM  DETECTION 
MALDEN  MA  (42°26'34"N,  71°05’07"W) 


TABLE  IB 

SECONDARY  SONIC  BOOM  DETECTION 
MALDEN  MA  (42°  26  »  31"N,  71°05,07"W) 


Degrees  re  True  North 

Hasped  Time  from  Measured  Event  to  Radar  Track  Time  CPA  Hyannis  Ma . 


TABLE  1C 

SECONDARY  SONIC  BOOM  DETECTION 
MALDEN  MA  (42°26'34"N,  71°05’07"W) 


CPA  Hynr 


TABLE  IE 

SECONDARY  SONIC  BOOM  DETECTION 
MISCELLANEOUS  SITES 


Degrees  relative  to  True  North. 

^Elapsed  Time  from  measured  Event  to  Radar  Track  Time  CPA  Hyannis, 


TABLE  2.  SECONDARY  SONIC  BOOM  REPORTS 


Date 


Observation  Time 


Locat ion 


Apr: 

LI  23 

4:00 

p .  m . 

Apr: 

LI  24 

3:55 

p .  m . 

May 

11 

5:17 

p.m. 

May 

22 

8:20 

a  .m . 

5:15 

p.m. 

May 

23 

8:23 

a .  m . 

9:41 

a  .m. 

May 

24 

8:22 

a .  m . 

May 

27 

7:45 

a .  m . 

8:23 

a .  m . 

9:36 

a  .m. 

9:40 

a .  m . 

9:40 

a  .m. 

9:40 

a .  m . 

May 

28 

8:16 

a .  m . 

May 

29 

8:35 

a  .m. 

9:09 

a  .m. 

9:36 

a .  m . 

9:37 

a  .m . 

May 

30 

9:08 

a .  m . 

9:34 

a .  m . 

9  to 

10  a 

8:15 

a .  m . 

May 

31 

7  :  30 

a  .m . 

8:12 

a  .m. 

8:15 

a  .m . 

Thompson  CT 

Thompson  CT 

Thompson  CT 

Milton  MA 
Thompson  CT 

Wilmington  MA 
Wilmington  MA 

Natick  MA 

Natick  MA 
West  Newton  MA 
West  Newton  MA 
Natick  MA 
Walpole  MA 
Wilmington  MA 

West  Newton  MA 

West  Newton  MA 
Natick  MA 
Brookline  MA 
Way land  MA 

Walpole  MA 
Walpole  MA 

m.  Medfield  MA 

Wilmington  MA 

Natick  MA 
Brookline  MA 
Walpole  MA 
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TABLE  3.  MEASURED  AND  COMPUTED  ARRIVAL  TIMES  OF  SECONDARY  SONIC  BOOM  SIGNATUTES  AT 
MALDEN  MA  SITE 
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TABLE  4.  MEASURED  AND  COMPUTED  ELEVATION  ANGLES  OF  INCIDENT 
SECONDARY  SONIC  BOOM  ARRIVALS  AT  MALDEN  MA  SITE 


Plot 

ID 

Flight 

, 1  Date 

Measured  Elevation 
Angle 
(Degrees) 

Computed  Elevation 
Angle2 
(Degrees) 

a 

BA 

5/23 

19 

7 

b 

BA 

5/29 

11 

7 

c 

AF 

5/30 

5 

Np3 

d 

BA 

5/31 

12 

6 

e 

AF 

6/13 

17 

23 

f 

BA 

6/13 

10 

24 

g 

BA 

6/20 

16 

12 

h 

AF 

6/27 

6 

22 

i 

BA 

6/27 

29 

23 

j 

AF 

7/18 

11 

10 

k 

BA 

7/18 

29 

17 

1 

AF 

7/25 

NM4 

17 

m 

BA 

7/25 

NM 

7 

1 

BA  171 

and  AF  001 

flights 

2  Based  on  ray  path  of  Type  I  arrival 

3  No  secondary  boom  predicted 

4  Not  measured 

(See  Figure  50  for  data  plot.) 
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APPENDIX  A 


TSC  SONIC  BOOM  COMPUTER  PROGRAM 

A. 1  EVOLUTION  OF  THE  PROGRAM 

A  program  made  available  to  TSC  for  the  present  study  was  a 

modification  by  E.  B.  Wright,  dated  April  1978,  of  an  "Atmospheric 

Raytrace  Program"  from  Sandia  Laboratories,  Albuquerque  NM. 

References  cited  in  the  transmitted  source  program  listing  were  a 

1 2 

1967  Sandia  Laboratories  report  by  Dr.  R.  J.  Thompson  and  a  1972 
article13  in  the  "Journal  of  the  Acoustical  Society  of  America," 
also  by  Thompson. 

Both  the  program  listed  in  Thompson's  report  and  Wright's 
modification  were  evidently  intended  for  analyses  of  acoustic 
arrivals  from  compact  explosive  sources  in  the  atmosphere.  A 
supersonic  aircraft  is  a  considerably  different  type  of  source,  but 
the  same  propagation  principles  apply  to  both  cases. 

In  Progress  Report  2  of  the  present  study, 1^  some  results 
were  reported  that  were  computed  with  a  somewhat  tedious  applica¬ 
tion  of  the  Wright  version  of  the  program  and  with  relatively 
crude  approximations.  Since  that  time,  a  new  program  has  been 
written  that  builds  upon  certain  key  subroutines  in  the  Thompson- 
Wright  program  and  that  is  specifically  intended  for  the  interpre¬ 
tation  of  secondary  sonic  boom  arrivals  from  inbound  flights  to 
JFK. 

A. 2  INPUT  DATA 

The  input  data  to  the  TSC  version  of  the  program  consists  of 
meteorological  data  and  flight  track  data.  The  format  in  which 
such  data  was  available  to  us  determined  the  choice  of  input  format 
for  the  program. 

Meteorological  Data.  The  Air  Resources  Laboratory  (Silver 
Springs  MD) ,  of  the  National  Oceanic  and  Atmospheric  Administra¬ 
tion  (NOAA)  furnished  TSC  with  synthesized  meteorological  soundings 


for  the  Boston  area  for  the  months  of  May,  June,  and  July.  The 
sources  and  procedures  for  the  synthesized  profiles,  as  explained 
in  a  letter  of  June  7,  1979  from  Albion  D.  Taylor,  are  as  follows. 

The  Upper  Air  Branch  archive  tape  records,  which  constituted 
one  source,  gives  heights  and  temperatures  for  each  point  of  a 
grid  (grid  spacing  approximately  381  km  or  206  nautical  miles)  cov¬ 
ering  the  Northern  Hemisphere  at  each  of  a  selected  set  of  pressure 
levels.  The  levels  ordinarily  used  are  those  of  1000,  850,  700, 
500,  400,  300,  250,  150,  and  100  mb.  The  corresponding  grid  values 
of  temperatures  and  heights  for  such  levels  are  derived  from  the 
fitting  of  smooth  fields  through  all  available  rawinsonde  (balloon 
borne)  observations;  this  fitting  is  customarily  carried  out  by 
the  National  Weather  Service  (NWS)  as  part  of  its  day-to-day 
operations.  Grid  values  for  levels  of  70,  50,  30,  and  10  mb  were 
obtained  in  nearly  the  same  way,  but  had  to  be  specially  analyzed 
by  the  Upper  Air  Branch. 

Data  corresponding  to  pressure  levels  of  5,  2,  1,  and  0.4  mb 
were  obtained  by  an  analysis  of  upper  atmospheric  temperatures 
as  deduced  from  the  remote  sensing  of  infrared  radiation.  The 
sensing  of  the  infrared  radiation  from  the  atmosphere  was  done  by 
instrumentation  carried  on  satellites. 

The  gridded  data  with  suitable  interpolations  yields  height 
and  temperatures  for  each  considered  pressure  level  for  the 
atmosphere  above  Boston  (whose  location  was  taken  as  71°W,  42.3°N). 
It  also  yielded  the  corresponding  values  (with  finite  difference 
computations)  for  the  two  horizontal  components  of  the  pressure 
gradient  for  each  pressure  level.  The  horizontal  pressure 
gradient  in  turn,  via  the  "geostrophic  wind  assumption",  yields 
the  horizontal  wind  components  at  that  same  point. 

The  computer  listing  supplied  to  us  by  NOAA  gives,  for  each 
day  and  for  1200  GMT  (0800  EDST) ,  the  following  quantities: 

a)  altitude,  Y(I),  in  meters  above  mean  sea  level 

b)  temperature,  T(I),  in  °C 
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c)  wind  speed  magnitude,  W(I),  in  m/s 

d)  direction,  D(I),  in  degrees,  reckoned  clockwise  from  true 
north,  from  which  the  wind  is  blowing. 

Here  the  index  I,  ranging  typically  from  1  to  18,  corresponds 
to  increasing  altitude,  Y(I  +  1)  >  Y(I),  and  to  decreasing  pres¬ 
sure.  The  symbols  used  here  are  those  that  appear  in  the  FORTRAN 
listing  of  the  computer  program. 

Flight  Track  Data.  Tabulated  data  obtained  by  the  air  traffic 
control  radar  system  for  inbound  Concorde  flights  into  JFK  were 
supplied  to  TSC  by  the  FAA  Data  Systems  Office  of  Nashua  NH.  These 
data  include  aircraft  flight  number,  speed  in  knots,  altitude, 
position,  and  time  of  day.  The  numbers  that  the  program  ordinarily 
requires  for  input  are 

a)  height,  HS(I),  in  hundreds  of  feet 

b)  X-coordinate ,  XS(I) ,  in  nautical  miles 

c)  Y-coordinate ,  YS(I),  in  nautical  miles 

d)  aircraft  ground  speed,  SMACH(I),  in  knots 

e)  time  of  day,  TAU(I),  generally  expressed  as  hour  (first 
two  digits),  minutes  (next  two  digits),  and  seconds  (next 
two  digits,  decimal  point,  digit).  The  FAA  tabulation 
gives  time  GMT  but,  for  the  summer  months,  we  ordinarily 
subtract  04  from  the  first  two  digits  to  obtain  time 
EDST. 

The  program  does  not  actually  use  the  input  ground  speed  but 
computes  it  internally  from  the  X,  Y  position  versus  time.  How¬ 
ever,  unless  at  least  one  ground  speed  larger  than  100  is  input, 
the  program  interprets  input  data  differently  than  as  discussed  in 
this  summary. 

FAA  personnel  in  Nashua  say  that  the  X  and  Y  coordinates  in 
the  FAA  printout  correspond  to  a  good  approximation  to  distance 
increments  in  nautical  miles  in  eastward  and  northward  directions, 
respectively.  The  origin  is  such  that  a  point  in  Hyannis  with 
latitude  and  longitude  41.7125°N,  70.21S6°W  is  at  X  =  365.6250  and 
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Y  =  218.3750.  The  conversion  suggested  by  these  numbers  is 
LON  =  70.2156  +  (365.6250  -  X)/44.790 
LAT  =  41.7125  +  (Y  -  2 18 . 3750 ) /60 . 000 

such  that  1  degree  longitude  corresponds  to  AX  =  -  44.790  nautical 
miles,  1  degree  latitude  corresponds  to  AY  =  60.000  nautical  miles. 
The  coordinate  origin  would  then  be  78.38  W,  38.07  N,  which  is  in 
the  vicinity  of  Charlottesville  VA. 

Input  Format .  A  typical  input  set  includes  the  following 
(see  Figure  A-l > : 

Line  1:  month  (integer) 

day  of  month  (integer) 
vear  (integer) 

1  if  airplane  is  flying  into  JFK 
1  if  airplane  is  British  Airways,  2  if  airplane 
is  Air  France 

1  if  airplane  is  BA,  2  if  airplane  is  Air  France 
flight  number  (integer) 

1  if  NOAA-synthes i zed  Boston  meteorological  data  is 
used. 

Line  2:  integer  NZ  giving  number  of  weather  data  altitudes 

NZ  Lines:  height,  temperature,  wind  direction,  wind  speed 
(lines  in  order  of  increasing  height) 

Next  Line:  integer  NS  giving  number  of  aircraft  positions 
along  flight  track 

NS  Lines:  aircraft  height,  X-coordinate ,  Y-coordinate ,  0.0, 
ground  speed,  and  time  of  day.  (The  0.0  is  input 
if  aircraft  headings  are  to  be  computed  from  flight 
trajectory.  Lines  must  be  in  order  of  increasing 
time. ) 

A. 3  REDUCTION  AND  SMOOTHING  OF  INPUT  DATA 

The  computations  produced  on  a  given  input  data  set  consist 
of  (1)  reduction  and  initial  processing  of  the  input  data,  (2) 
prediction  of  secondary  boom  impact  areas,  (3)  prediction  of  event 


arrival  times  at  a  specified  reception  site,  and  (4)  prediction  of 
extent  of  primary  carpet  region.  The  first  such  category  is  dis¬ 
cussed  below. 

Reduced  Meteorological  Data.  The  initial  steps  in  the  program 
transform  the  input  meteorological  data  to  a  tabulation  that  is 
generally  longer  (35  lines)  and  smoother  than  the  probable  atmos¬ 
pheric  profile  over  Boston.  Such  a  tabulation  appears  in  the  output 
(Figure  A-2)  and  gives  altitude  in  meters,  height  in  meters,  wind 
speed  in  meters  per  second,  direction  from  which  wind  is  blowing  in 
degrees,  and  the  sound  speed  in  meters  per  second.  The  height  is 
simply  altitude,  adjusted  such  that  the  initial  height  is  zero.  The 
sound  speed  is  computed  from  the  input  temperature  (after  smoothing) 
using  the  standard  theoretical  relation  for  the  sound  speed  in  air. 

The  reasons  for  the  smoothing  of  the  input  data  are  (1)  the 
accuracy  of  the  NOAA  predictions  for  the  time  of  the  event  is 
probably  far  less  than  any  roughness  implied  by  the  input;  (2) 
details  of  such  roughness  even  if  real  would  not  be  expected  to  be 
the  same  all  along  the  acoustic  propagation  path;  (3)  non-physical 
roughness  in  the  meteorological  profile  will  lead  to  spurious 
predictions  of  acoustical  propagation  that  would  confuse  the 
interpretation  of  the  calculations. 

Reduced  Flight  Trajectory  Data.  The  input  X(t)  and  Y(t)  tab¬ 
ulations  are  smoothed  to  eliminate  random  errors  in  the  radar  track¬ 
ing.  The  ground  speed  vector  is  then  computed  by  numerical  differ¬ 
entiation,  yielding  dX/dt  and  dY/dt  for  each  of  the  considered  times 
along  the  flight  trajectory.  Subtraction  of  the  wind  velocity  at  the 
corresponding  flight  altitude  produces  the  airplane's  velocity  rel¬ 
ative  to  that  of  the  air.  The  direction  of  that  vector  gives  the 
aircraft  heading  HDG(I);  division  of  its  magnitude  (air-speed)  by 
the  flight  altitude  sound  speed  gives  the  Mach  number ,  which  there¬ 
after  is  referred  to  by  the  program  as  SMACH(I).  The  X  and  Y  coor¬ 
dinates  are  converted  to  longitude  and  latitude  using  the  equation 
given  in  Section  A. 2. 

The  present  version  of  the  program  is  written  for  inbound 
flights  and  computes  the  aircraft  heading,  assuming  it  is  between 
180°  (southward)  and  270°  (westward).  Also,  the  program  assumes 


that  the  initial  flight  track  point  corresponds  to  supersonic 
flight  and  disregards  all  flight  track  data  after  the  speed  first 
goes  subsonic  (if  ever). 

Another  weeding-out  process  eliminates  the  first  point  (and 
then  all  subsequent  ones)  for  which  the  Mach  number  and  aircraft 
heading  are  such  that  none  of  the  sound  rays  generated  toward  the 
right  of  the  airplane  reach  the  ground.  The  theoretical  model  is 
such  that  a  ray  leaving  the  airplane  with  a  bearing  angle  0g  in 
the  horizontal  direction  must  have  an  initial  elevation  angle 
consistent  with  the  local  geometry  of  the  Mach  cone  relative  to 
the  wind.  This  requires  in  particular  that  the  horizontal  phase 
velocity  Vp  of  the  ray  is  =  u^  +  Mc^  cos  (8g  -  ©H) ,  where  u^  is 
the  wind  velocity  component  (positive  if  blowing  towards  listener) 
in  direction  0g  at  the  altitude  of  flight,  c^  is  the  sound  speed, 
0^  is  the  aircraft  heading,  and  M  is  the  Mach  number. 

The  acoustical  model  predicts  that  Vp  is  constant  along  any 

given  ray  path.  A  ray  with  horizontal  phase  velocity  can  reach 

the  ground  only  if  V  >  u  +  c  ,  where  u  and  c  are  the  wind 
s  7  P  g  g  g  g 

velocity  component  (direction  0g)  and  sound  speed  at  the  ground. 

A  flight  trajectory  data  point  is  considered  a  potential  source  of 

ground  level  sonic  booms  only  if  this  inequality  is  satisfied  for 

at  least  one  bearing  angle  0g. 

With  the  application  of  the  criteria  just  described,  the 
number  of  source  points  NS  is  in  general  reduced.  The  program 
prints  out  a  listing  (Figure  A-3)  of  aircraft  height  in  meters, 
longitude  and  latitude  of  airplane,  heading  angle  of  airplane,  and 
Mach  number,  versus  increasing  time  (hours,  minutes,  seconds 
format)  for  each  input  flight  trajectory  point  that  survives  this 
scrutiny. 

A. 4  SECONDARY  SONIC  BOOM  FOOTPRINTS 

From  each  source  point  along  the  aircraft  trajectory,  the 
program  traces  up  to  21  rays  that  carry  secondary  sonic  booms  from 
the  aircraft  to  the  ground.  Two  types  of  rays  are  distinguished: 
(1)  rays  that  go  from  the  airplane  upwards  to  mesosphere  and  then 
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back  to  the  ground  and  (2)  rays  that  go  from  the  airplane  downwards 
to  ground  (as  a  primary  boom),  that  are  then  ground-reflected, 
propagate  up  to  the  mesosphere,  then  return  to  the  ground. 

Neither  type  of  ray  is  possible  unless  the  profile  of  sound 
speed  and  wind  succeeds  in  refracting  the  ray  downward  in  the 
atmosphere.  The  criterion  for  this  to  occur  for  a  given  ray  bear¬ 
ing  angle  0g  is  that 


V_ 


(u  +  c) 


max 


where  (u  +  c)max  is  the  maximum  value  of  sum,  wind  velocity  com¬ 
ponent  in  direction  0g  plus  sound  speed,  that  occurs  among  heights 
exceeding  the  height  of  the  aircraft.  Here  Vp  is  the  phase 
velocity  defined  in  the  preceding  subsection.  Combining  the 
relations  given  there  with  that  above  yields. 


ug  +  cg  <  Uf  +  Mc£  cos  (eB  -  eH)  <  (u  ♦  c)max 

It  should  be  noted  that  ug,  uf  and  (u  +  c)max  depend  on  the  bearing 
angle. 

In  general,  there  is  only  a  narrow  range  (if  any  at  all)  of 
bearing  angles  for  which  the  above  criteria  are  satisfied,  this 
range  being  typically  of  the  order  of  5°.  What  the  program  does 
is  to  trace  out  up  to  21  rays  from  the  aircraft  that  have  bearing 
angles  within  this  range.  For  each  ray  of  a  given  type  (I  or  II 
as  defined  above),  the  computation  yields  a  horizontal  distance 
over  which  the  ray  propagates  from  aircraft  to  ground  and  a  time 
of  transit.  Results  for  each  type  of  ray  and  for  each  source 
location  are  printed  out  in  the  manner  indicated  in  Figure  A-4. 

The  title  on  the  printout  indicates  the  considered  secondary 
ray  type,  (I)  airplane  to  stratosphere  (mesosphere)  to  ground  or 
(II)  airplane  to  ground  to  stratosphere  (mesosphere)  to  ground. 

Then  source  location  (airplane  altitude,  longitude,  and  latitude) 
is  identified,  along  with  the  airplane  heading,  Mach  number,  and 
the  time  (hours,  minutes,  and  seconds)  at  which  the  airplane 


passed  through  this  location.  Two  other  parameters  indicate 
where  the  Malden  monitoring  site  is  with  respect  to  this  source 
location.  The  bearing  of  the  vector  pointing  from  source  to 
Malden  is  in  degree  clockwise  from  true  north.  The  distance  to 
Malden  is  printed  in  kilometers,  where  1  km  =  0.54  nautical  miles. 

In  the  accompanying  ray  tabulation,  the  computed  rays  are 
numbered  in  the  order  of  increasing  bearing  angle.  For  each  such 
ray,  the  ray  range  (horizontal  distance  travelled)  is  given  in 
kilometers.  This  range  typically  begins  large,  decreases  to  a 
minimum,  then  increases  to  a  slightly  higher  value,  when  the 
bearing  angle  runs  through  the  range  of  possible  values.  The 
reason  for  this  behavior  is  well  understood  in  terms  of  the  over¬ 
all  ray  paths.  The  minimum  coincides  with  where  the  ray  caustic 
(surface  where  adjacent  rays  intersect)  touches  the  ground. 

The  tabulation  also  gives  the  longitude  and  latitude  of  the 
point  at  which  each  ray  touches  the  ground  and  the  time  (hours, 
minutes,  and  seconds)  at  which  this  occurs.  (Such  computations 
of  course  depend  on  the  location  of  the  source  and  the  time  of  ray 
launching.)  The  transit  time  of  the  ray  along  the  path  is  in  a 
minutes  and  seconds  format.  The  miss  distance  printed  is  the 
distance  of  the  ray  arrival  point  from  the  Malden  site. 

The  line  traced  out  on  the  earth's  surface  by  the  ray  arrival 
points  from  any  one  such  tabulation  gives  what  is  here  termed  a 
"footprint."  The  airplane  generates  one  such  footprint  of  each 
type  from  each  point  along  its  flight  trajectory.  The  set  of  all 
footprints  of  a  given  ray  type  (I  or  II)  covers  an  area  of  the  map 
which  roughly  coincides  with  where  the  secondary  booms  of  that 
type  are  expected  to  impact.  The  density  of  rays  (rays  per  unit 
area  of  ground  surface)  gives  a  rough  indication  of  how  intense 
the  arrivals  may  be  expected  to  be.  Since  the  density  is  higher 
for  portions  of  the  impacted  region  that  are  closer  to  the  flight 
track,  the  quantities  of  greatest  interest  are  the  focus  lines. 

The  focus  line  for  type  I  rays  is  always  closer  to  the  flight  track 
than  that  for  type  II  rays.  The  program  does  not  actually  deter¬ 
mine  and  plot  these  inner  envelope  lines,  but  such  are  easily  de¬ 
rived  by  hand-plotting  a  sequence  of  ray  footprints  on  a  map. 


A. 5  RECEPTION  TIME  COMPUTATIONS 


For  each  secondary  boom  ray  type,  the  program  endeavors  to 
find  the  ray  that  comes  closest  to  the  Malden  site.  All  the  rays 
successively  computed  are  compared  according  to  their  "miss  dis¬ 
tances"  and  detailed  intermediate  results  are  kept  stored  for  the 
closest  ray  found.  Successive  iterations  give  progressively 
closer  rays. 

Figure  A-5  shows  typical  output  regarding  such  a  "closest" 
ray.  The  listing  gives  that  ray's  miss  distance  and  detailed 
information  concerning  the  point  on  the  aircraft  trajectory  from 
which  the  ray  was  launched.  Also  listed  is  the  closest  ray's 
range  (horizontal  distance  traveled) ,  transit  time  and  bearing 
angle. 

The  interpretation  of  such  output  depends  on  whether  the 
Malden  site  falls  on  the  near  side  or  far  side  of  the  corresponding 
inner  envelope  line.  If  it  falls  on  the  far  side  (typical  for  the 
type  I  secondary  boom  carpet),  then  further  narrowing  of  the  inter¬ 
vals  between  successively  considered  points  along  the  aircraft 
trajectory  and  an  interpolation  will  in  principle  yield  a  ray 
(possibly  more  than  one)  that  connects  the  trajectory  with  the 
Malden  site.  In  such  a  case,  the  information  of  most  interest 
include  (1)  the  time  the  ray  arrives  at  Malden,  (2)  the  direction 
the  ray  appears  to  be  heading  when  it  reaches  Malden,  and  (3)  the 
elevation  angle  which  the  incoming  ray  makes  with  the  ground. 

These  are  parameters  that  can  be  compared  with  numbers  extracted 
from  analysis  of  the  received  signals  at  Malden.  Secondary  quan¬ 
tities  computed  by  the  program  (Figure  A-6)  are  the  ray’s  phase 
velocity,  the  angle  it  makes  with  the  horizontal  when  it  is 
launched,  and  the  maximum  height  the  ray  reaches  before  it  turns 
back  towards  the  ground.  The  program  also  lists  the  complete 
trajectory  for  that  ray;  range  and  height  versus  time.  Also 
listed  is  the  effective  atmospheric  profile  for  the  corresponding 
bearing  angle  (Figure  A- 7 ) :  U(Z)  is  the  wind  speed  component 
blowing  parallel  to  vertical  plane  of  ray  towards  Malden,  S(Z)  is 
the  effective  sound  speed  C(Z)  +  U(Z). 
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On  the  other  hand,  if  the  Malden  site  lies  on  the  near  side  of 
the  inner  envelope  line  (typical  for  the  type  II  secondary  boom 
carpet),  the  smallest  miss  distance  is  an  estimate  of  how  far  the 
site  is  from  that  line.  The  refined  theory,  in  which  geometrical 
acoustics  is  modified  to  include  full-wave  effects,  predicts 
arrivals  on  the  near  side  of  the  envelope  line  as  well  as  on  the 
far  side,  but  predicts  that  such  arrivals  are  weaker  the  further 
away  from  the  line  the  listener  falls.  Consequently,  the  computed 
minimum  miss  distance  should  be  negatively  correlated  with 
amplitudes . 

An  estimated  time  of  arrival  of  type  II  secondary  booms  at 
Malden  when  the  site  lies  between  the  two  envelope  lines  (the 
usual  case)  is  computed  with  the  assumptions  that:  (1)  the  type 
II  arrival  has  the  same  launch  point  and  azimuth  as  does  the  type 
I  arrival;  (2)  the  time  difference  between  type  I  and  type  II 
arrival  times  is  independent  of  range  and  may  therefore  be  ex¬ 
tracted  from  computed  ray  arrival  times  at  a  point  further  on 
beyond  Malden  where  both  rays  are  actually  received. 

A. 6  PRIMARY  CARPET 

The  primary  carpet  generated  by  a  supersonic  flight  consists 
of  all  points  that  receive  direct  rays  from  the  aircraft  trajectory 
without  intervening  excursions  to  the  mesosphere.  The  rays  gener¬ 
ated  at  a  given  point  on  the  flight  trajectory  that  are  proceeding 
obliquely  downwards  reach  the  ground  only  if 

u(z)  +  c(z)  <  uf  +  Mcf  cos(6g  -  eH) 

for  all  heights  z  that  are  below  the  height  of  the  airplane.  Here 
the  various  symbols  are  as  defined  in  Section  A. 3;  the  right  side 
of  the  inequality  is  also  the  phase  velocity  V  . 

A  ray  satisfying  the  criteria: 

Ug  +  Cg  =  Uf  +  MCf  C°S  (0B  '  eH3 
u ( z )  +  c  (z)  <  Ug  +  cg;  0  <z<  zf 

(d/dz)  (u(z)  +  c(z))  0;  z  =  0 
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is  such  that  it  propagates  to  the  ground  and  there  just  grazes  the 
ground.  Its  phase  velocity  marks  the  borderline  between  where  rays 
reflect  and  where  they  bend  back  upwards  before  they  touch  the 
ground.  The  point  at  which  it  touches  the  ground  should  therefore 
lie  on  the  border  between  where  rays  touch  the  ground  and  where 
they  do  not  touch  the  ground;  the  former  region  is  the  primary 
carpet. 

In  the  program’s  computation  of  the  primary  carpet  border 
(see  Figure  A-8),  the  first  of  the  three  equations  above  is  solved 
numerically  for  0g.  The  value  of  found  (should  more  than  one 
exist)  is  that  corresponding  to  the  smallest  positive  0g  -  0^  (or 
the  northwest  side  of  the  primary  carpet  for  inbound  flights). 

The  present  version  of  the  program  does  not  test  whether  the  other 
two  criteria  are  met;  the  input  meteorological  data  is  such  that 
the  check  need  not  be  made. 

Once  0g  is  found,  the  program  calculates  the  corresponding 
limiting  ray,  from  aircraft  trajectory  to  carpet  edge.  The  tabu¬ 
lation  (Figure  A-8)  lists,  for  each  aircraft  position,  the  ray 
launch  time  (hours,  minutes,  and  seconds),  which  is  the  time  the 
airplane  passes  through  the  point  from  which  the  ray  under  consid¬ 
eration  was  launched.  The  tabulation  gives  the  bearing  angle  of 
the  limiting  ray  and  the  distance  (km)  it  travels  before  it 
touches  the  ground.  The  longitude  and  latitude  of  that  point  is 
also  printed,  as  well  as  the  time  that  it  takes  the  ray  to  traverse 
that  distance. 

Other  columns  in  the  tabulation  relate  to  the  estimation  of 
the  properties  of  the  ground  wave  that  passes  the  Malden  site. 

The  ground  wave  is  a  very  weak  arrival  that  travels  to  the  edge  of 
the  primary  carpet  along  the  geometrical  acoustics  limiting  ray 
and  that  thereafter  travels  with  the  same  bearing  along  the 
ground  with  speed  c  +  u  Thus  each  limiting  ray  to  the  primary 

O  6 

carpet  edge  generates  a  ground  wave  ray  that  travels  along  the 
ground  in  the  same  general  direction.  One  (or  more)  such  ground 
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wave  ray  passes  through  the  Malden  site.  The  tabulation  in 
Figure  A-8  consequently  lists  the  distance  by  which  each  computed 
ground  wave  path  misses  Malden.  Interpolation  consequently  will 
give  the  time  at  which  the  ground  wave  passing  through  the  Malden 
site  was  launched. 

An  estimate  of  ground  wave  arrival  time  at  the  Malden  site 
comes  from  an  extended  version  of  Fermat's  principle  of  least  time. 
The  principle,  which  is  founded  on  basic  acoustical  theory,  states 
that,  if  one  takes  any  given  point  on  the  flight  trajectory  and 
computes  the  ficticious  arrival  time  FERMAT,  equal  to 

FERMAT  =  LNCHTM  +  RAYTM  +  EXTMN 

then  FERMAT  is  later  than  the  actual  ground  wave  arrival  time  at 
Malden.  Here  LNCHTM  is  the  time  the  grazing  ray  was  launched, 

RAYTM  is  the  transit  time  along  the  grazing  ray  to  the  primary 
carpet  edge,  and  EXTMN  is  the  time  it  would  take  sound  to  travel 
in  a  straight  line  (skewed  with  respect  to  the  grazing  ray)  along 
the  ground  from  the  grazing  ray  touch-down  point  to  the  Malden 
site.  The  speed  along  the  latter  leg  of  the  hypothetical  path  is 
that  appropriate  to  the  bearing  angle  from  the  touch-down  point 
to  Malden. 

The  minimum  value  of  the  computed  value  of  FERMAT  in  the 
tabulation  yields  an  estimate  of  the  actual  ground  wave  arrival 
time  in  Malden  and  also  corresponds  to  the  smallest  miss  distance. 
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FIGURE  A- 2.  REDUCED  METEOROLOGICAL  DATA  PRINTOUT 
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FIGURE  A- 3.  REDUCED  FLIGHT  TRAJECTORY  DATA  PRINTOUT 
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FIGURE  A- 4.  SAMPLE  SECONDARY 
DETERMINATION  OF  TYPE  I  FOCUS 
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FIGURE  A- 5.  SAMPLE  OUTPUT  REGARDING  AIRCRAFT  POSITION  AND  FLIGHT 
PARAMETERS  AT  TIME  OF  LAUNCHING  OF  RAY  THAT  GOES  TO  MEASUREMENT 
SITE 
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FIGURE  A-6.  SAMPLE  OUTPUT  REGARDING  RAY  TOUCHING  GROUND 
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FIGURE  A- 7.  SAMPLE  OUTPUT  GIVING  EFFECTIVE  SOUND  SPEED  PROFILE 
FOR  PROPAGATION  OF  RAYS  FROM  FLIGHT  TRACK  TO  MEASUREMENT  SITE 
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FIGURE  A- 8.  SAMPLE  OUTPUT  DESCRIBING  BORDER  OF  SONIC  BOOM 
PRIMARY  CARPET 
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